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ABSTRACT 

We combine near-UV spectra obtained with the Hubble Space Telescope GHRS 
echelle with far-UV spectra obtained with IMAPS and Copernicus to study the abun- 
dances and physical conditions in the predominantly ionized gas seen at high velocity 
(—105 km s~^ ^ Vq < —65 km s~^) and at intermediate velocity (—60 km ^ Vq < 
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—10 km s"-*^) along the line of sight to the star ( Ori. We have high resolution (FWHM 
~ 3.3-4.5 km s~"^) and/or high S/N spectra for at least two significant ions of C, N, 
Al, Si, S, and Fe — enabling accurate estimates for both the total A^(H II) and the 
elemental depletions. C, N, and S have essentially solar relative abundances; Al, Si, 
and Fe appear to be depleted by about 0.8, 0.3-0.4, and 0.95 dex, respectively, relative 
to C, N, and S. While various ion ratios would be consistent with collisional ionization 
equilibrium (CIE) at temperatures of 25,000 80.000 K, the widths of individual high- 
velocity absorption components indicate that T r-^ 9000ib2000 K — so the gas is not 
in CIE. Analysis of the C II fine-structure excitation equilibrium, at that temperature, 
yields estimates for the densities (rig ~ nn ~ 0.1-0.2 cm^^), thermal pressures {2n}iT 
~ 2000-4000 cm~^K), and thicknesses (0.5-2.7 pc) characterizing the individual clouds. 
We compare the abundances and physical properties derived for these clouds with those 
found for gas at similar velocities toward 23 Ori and r CMa, and also with several dif- 
ferent models for shocked gas. While the shock models can reproduce some features of 
the observed line profiles and sonic of the observed ion ratios, there are also significant 
differences between the models and the data. The measured depletions suggest that 
roughly 10% of the Al, Si, and Fe originally locked in dust in the pre-shock medium 
may have been returned to the gas phase, consistent with recent predictions for the 
destruction of silicate dust in a 100 km s~^ shock. The observed near-solar gas phase 
abundance of carbon, however, appears to be inconsistent with the predicted longer 
time scales for the destruction of graphite grains. 

Subject headings: ISM: abundances — ISM: general — stars: individual {( Ori, r CMa) 
— ultraviolet: ISM 

1. Introduction 

Interstellar gas found at high and intermediate velocities in the Galactic plane near associations 
of young, massive stars presumably reflects the combined effects of stellar winds and supernovae 
arising in those associations. Observations of UV absorption lines can reveal the kinematics, elemen- 
tal abundances (and depletions), ionization state, and physical properties (density, temperature, 
etc.) characterizing such high- and intermediate- velocity gas — and can thus provide significant 
empirical tests for models of (1) the interaction of supernovae with the surrounding interstellar 
medium, (2) shocks and the non-equilibrium evolution of shocked gas, and (3) the survival of dust 
grains in such environments. 

UV spectra obtained with Copernicus revealed absorption from various singly and doubly ion- 
ized species at both high (—100 to —80 km s~^) and intermediate (—80 to —15 km s~^) heliocentric 
velocities toward a number of stars in Orion (Cohn & York 1977; Cowie, Songaila, & York 1979). 
Cowie et al. ascribed the widespread high-velocity (HV) gas to a radiative shock produced by 
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one or more recent supernovae and called it "Orion's Cloak". Subsequent UV spectra of addi- 
tional stars in Orion obtained with lUE provided further information on the distribution of the 
HV and intermediate-velocity (IV) gas throughout that region (Shore 1982; Sonneborn, Shore, & 
Brown 1988). More recently, higher resolution UV spectra of 23 Ori obtained with the Hubble 
Space Telescope Goddard High Resolution Spectrograph (HST/GURS) have enabled more precise 
characterization of the HV and IV gas in that line of sight (Trapero et al. 1996; Welty et al. 1999). 
Toward 23 Ori, the HV gas is essentially fully ionized, and is characterized by a density Ue ~ nu 
~ 0.4-0.5 cm^^, a temperature T ~ 8000 K, and a thermal pressure 2n}iT 5-10 x 10^ cm~^K. 
As the ion ratios A^(X III)/A^(X II) for aluminum and silicon would be consistent with collisional 
ionization equilibrium (CIE) at much higher temperatures — of order 25000 K — the HV gas is 
clearly not in CIE. No data (or only upper limits) were available, however, for a number of poten- 
tially significant ions (e.g., C III, N III, S III, Fe II, Fe III), and some of the spectra were obtained 
at only moderate spectral resolution. 

In this paper, we discuss the predominantly ionized HV and IV gas toward a second Orion 
star {( Ori), using near-UV spectra obtained with GHRS and far-UV spectra obtained with the 
Interstellar Medium Absorption Profile Spectrograph (IMAPS) and Copernicus. In Sections 2 and 3, 
we describe the spectra and our analysis of the observed absorption-line profiles. In Section 4, we 
discuss the abundances and physical properties derived for the individual components in the HV 
and IV gas. The unprecedented combination of high spectral resolution (FWHM ~ 3.3-4.5 km 
s"-*^), high signal-to-noise ratio (S/N), and broad spectral coverage has enabled the determination 
of accurate abundances for both singly and doubly ionized C, N, Al, Si, S, and Fe — which 
then provide strong constraints on the overall A^(H), the depictions, and the ionization in the HV 
and IV gas. Detailed multicomponent fits to the high-resolution line profiles and analyses of the 
C II fine-structure excitation equilibrium yield values for the temperature and density of the gas, 
respectively. In Section 5, we compare the elemental abundances and derived physical properties 
with those found for gas observed at similar velocities in several other lines of sight, we compare the 
observed line profiles and ion ratios [A'"(X III)/A^(X II), N{X IV)/N{X II)] with those predicted by 
various models for ionized gas (including newly-computed models for shocked gas), and we discuss 
the observed depletions in relation to recent models for dust processing in shocks. In an appendix, 
we present similar (though more limited) data for the high-velocity gas toward r CMa. 



2. Observations and Data Reduction 

2.1. C Ori Line of Sight 

The 09.5 lb star C Ori A [V = 1.88; E{B-V) = 0.06; v,^d = 18.1 km s'^; vsini = 135 km s"^] 
is located at (/, b) = (206?45, — 16?59), at the southeast end of the belt of Orion. (Unless otherwise 
specified, we will use heliocentric velocities; the zero point for velocities with respect to the Local 
Standard of Rest is at v© ~ -1-17.6 km s~^.) The parallax 3.99±0.79 mas measured by Hipparcos 
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implies a distance of 251I42 pc (ESA 1997), consistent with, but perhaps somewhat smaller than 
the value 340ib85 pc inferred from the spectral type and apparent magnitude (using the absolute 
magnitudes of Blaauw 1963). Brown, Walter, &: Blaauw (1998) have suggested that Orion OB 
association subgroup lb (of which Ori is a member) lies at an average distance of about 430 pc, 
however, based on Hipparcos parallaxes for a number of Orion stars. The line of sight to ^ Ori 
has a total neutral hydrogen column density log[iV(H I + 2H2)] = 20.43 cm~^ (Bohlin, Savage, 
& Drake 1978; Shull & van Stecnbcrg 1985; Diplas & Savage 1994), with only a small fraction in 
molecular form (log[A^(H2)] 15.85 cm~^; Jenkins & Peimbert 1996). C Ori has been observed in 
connection with a number of surveys of optical and UV interstellar absorption- lines (e.g., Hobbs 
1969; MarschaU & Hobbs 1972; van Steenberg & Shuh 1988; Welty, Hobbs, k Kulkarni 1994; Welty, 
Morton, & Hobbs 1996), some properties of the primarily neutral gas at low velocities have been 
discussed by Drake & Pottasch (1977) and by Jenkins &; Peimbert (1996), and the high-velocity 
ionized gas associated with Orion's Cloak has been described by Cowie et al. (1979). 



2.2. GHRS Spectra 

We observed C Ori with the HST GHRS during Cycle 5 (on 1996 Feb 12), as part of GO 
program 5896 for studying high-velocity ionized gas. All observations reported here were obtained 
through the small science aperture, in order to achieve the highest possible spectral resolution. 
Eight wavelength regions between 1188 and 1676 A and six regions between 1740 and 2805 A, 
each covering 6-15 A, were observed using the echelle grating at resolutions of 3.3-3.7 km s 
(Table 1). In order to obtain proper sampling and to reduce the effects of photocathode sensitivity 
variations, all exposures were divided into a number of slightly shifted sub-exposures, using FP- 
SPLIT = DSFOUR and COMB = 4 (see, e.g., Duncan k Ebbets 1989). Approximately 6% of each 
exposure was devoted to measuring the inter-order background. Because the cores of many of the 
obviously saturated absorption lines at low velocities fell either a few percent above or below the 
zero level in the extracted spectra received from STScI, the spectra were re-extracted using the 
STSDAS CALHRS routines, with a cubic polynomial fit to the background (instead of the default 
constant background). Even after reprocessing, the cores of some saturated lines still differed 
slightly (generally less than 1%) from zero; those spectra were adjusted accordingly. In addition, 
some G160M spectra of ( Ori (originally acquired under GO program 2584) were obtained from 
the HST archive in order to measure absorption lines from species not covered in our spectra. Nine 
wavelength regions between 1186 and 1875 A, each covering approximately 36-38 A, had been 
observed with grating G160M at intermediate resolutions of 13-20 km (Table 1). 

Because ( Ori is very bright, even these short GHRS exposures (3.5-5.5 min) were expected 
to yield S/Ns of 100-500 per half resolution element in the extracted spectra — if the detector 

fixed-pattern noise (FPN) could be adequately removed. For the ECH-A and ECH-B spectra, the 
STSDAS routines POFFSETS, DOPOFF, and SPECALIGN were initially used to determine the 
detector fixed-pattern features and to align and sum the sub-exposures. Near most of the absorption 
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features in the spectra, application of those routines did yield S /N nearly as high as those expected 
from photon statistics. Near the extreme ends of some of the spectra, and in cases where the 
fixed-pattern features differed in strength among the sub-exposures, however, those routines did 
not completely remove the FPN. An alternative FPN removal procedure based on Fourier transform 
methods (developed by EBJ) was therefore also used to correct and combine the sub-exposures. 
In general, the two procedures yielded very similar combined, corrected spectra. The STSDAS 
routines performed better near the cores of the strongest absorption features, while the Fourier 
technique was more successful in removing FPN features which differed in strength among the 
sub-exposures (e.g., near the HV absorption for the Fe II A2600 line). 

For the G160M exposures obtained from the archive, the relative offsets among the subexpo- 
sures in the nominal wavelength scales supplied by STScI were determined by fitting the profiles 
of interstellar and/or stellar lines in the spectra. The individual subexposures were shifted accord- 
ing to those relative offsets (which were always less than 13 km s~^), sine interpolated (Spitzer & 
Fitzpatrick 1993) onto a roughly twice-oversampled standard wavelength grid, and summed. Ob- 
vious deviant pixels in the individual subexposures (due to detector pattern features, bad pixels, 
radiation hits, etc.) were given zero weight in the sums; no other explicit detector pattern noise 
identification and removal was performed, however. Adjacent pixels were then summed to yield 
approximately optimal sampling (i.e., two points per resolution element) in the final spectra. 

Portions of the spectra near interstellar absorption lines of interest (usually including 150-^ 
250 points) were normalized by fitting the continuum regions with Legendre polynomials of orders 
typically 4-9. The continuum signal-to-noise ratios achieved for the echelle spectra range from 
100-550 (Table 1); the corresponding equivalent width uncertainties {2a) for weak, unresolved 
absorption features range from 0.2-1.0 mA. For the G160M data, the S/N are typically between 
60 and 350, yielding 2a equivalent width upper limits or uncertainties of 1-5 mA. For some lines 
with strong absorption over a wide velocity range (e.g., C II A1334, Si III A1206, and the stronger 
Si II lines), the adopted continuum is less well determined. The increased uncertainty in continuum 
placement affects primarily the absorption at low and intermediate velocities. The quoted equivalent 
width errors include contributions from both photon noise (Jenkins et al. 1973) and continuum 
placement (Sembach Sz Savage 1992), added in quadrature. 

Some of the normalized GHRS spectra, together with similar spectra from IMAPS (§ 2.3) and 
Copernicus (§ 2.4), are shown in Figure 1. This figure includes the absorption at both low and 
high velocities; the high- and intermediate-velocity absorption features (as defined in § 3) will be 
shown in more detail in Sections 5.3.1 and 5.3.2, where they are compared with those predicted 
by several shock models. The equivalent widths (or limits) for the high-velocity absorption due to 
various atomic and ionic species found in the UV spectra are listed in Table 2. The intermediate- 
velocity absorption features are often blended with the strong absorption due to primarily neutral 
gas at lower velocities. The source of each observed line (GHRS echelle, GHRS G160M, IMAPS, 
Copernicus) is noted following the wavelength. Dominant contributors to blended lines are marked 
with a "B" following the equivalent width; minor contributors to blended lines are included in the 
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table, but are referred to the dominant contributors for the total equivalent widths. 

2.3. IMAPS Spectra 

The far-UV spectrum of ( Ori between about 950 and 1150 A was observed with IMAPS during 

that instrument's first orbital flight on the ORFEUS-SPAS platform, carried by Space Shuttle 
flight STS-51 in 1993 September. Jenkins et al. (1996) have given a very detailed description of 
the IMAPS payload, its performance during that orbital flight, and the data reduction procedures 
employed to obtain wavelength calibrated, extracted spectra from the numerous two-dimensional 
partial echelle images typically acquired for each target. Sixty-three CCD frames, each covering 
about 1/4 the free spectral range of the echelle, were obtained in 2412 seconds of total integration on 
Ori. These images were corrected for electronic interference effects, dark current, non-uniformities 
in photocathode response, and geometric distortion as described by Jenkins et al., and the individual 
frames were aligned and co-added. The one-dimensional spectra were obtained via an optimal 
extraction routine which took into account the contamination from the wings of adjacent orders. 

The resolution achieved in the IMAPS spectra of Ori is approximately 4.5 km (FWHM; 
Jenkins & Peimbert 1996). We therefore smoothed the original oversampled spectra, by adding 
adjacent data points, to obtain (roughly) two pixels per resolution element. The continua near lines 
of interest were fitted with low-order Legendre polynomials. The S/Ns measured in the IMAPS 
spectra are generally between about 30 and 80 for wavelengths between 1000 and 1150 A, but can 
be significantly smaller below 1000 A. Background levels are typically accurate to within about 5% 
(Sofia & Jenkins 1998); corrections could be made in cases where there are nearby saturated low- 
velocity absorption features. The velocity zero point was established by using the various telluric 
O I* and O I** features present in the spectrum. 

Jenkins &; Peimbert (1996) have discussed the numerous H2 absorption features arising in the 
low-velocity gas toward ( Ori. In this study of the ionized gas, we focus on the absorption lines of 
C II, C III, N II, N III, S III, S IV, and Fe III found in the IMAPS spectra. Several of the normalized 
IMAPS spectra are shown in Figure 1 and the equivalent widths measured from those spectra are 
listed in Table 2. The quoted uncertainties include contributions from photon noise, continuum 
placement, and background level — added in quadrature. 

2.4. Copernicus Spectra 

Relatively high S /N far-UV spectra of ^ Ori had also been obtained with Copernicus by Cowie 
et al. (1979) for their study of the ionized high- velocity gas in the Orion region. Although of 
somewhat lower resolution (FWHM ~ 0.045 A, corresponding to 13.5 km s~^ at 1000 A), these 
spectra were retrieved from the Copernicus archive^ as a supplement to the IMAPS spectra in 
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cases where the latter were of low S/N (e.g., for C III and N III). The Copernicus spectra were 
also useful for confirming the weak absorption due to S IV (at low velocities) and Fe III (at high 
velocities) in the IMAPS spectra. In general, we chose the group of contemporaneous UI scans 
with the widest spectral coverage encompassing each line of interest — to minimize thermal shifts 
and sensitivity variations and to reach the continuum on either side of the absorption feature. 
The individual scans were co-added using a version of the Copernicus STACK program; continua 
were fitted using low-order Legendre polynomials. Background levels could be established for most 
lines either via reference to the typically saturated absorption at low velocities or by comparison 
with the corresponding IMAPS spectra. The normalized spectra of C III and N III are shown in 
Figure la (with source name OAO) and the equivalent widths measured from those spectra are 
listed in Table 2. These spectra exhibit very good agreement with those shown by Cowie et al. 
(1979) — which suggests that we have used the same sets of scans for each line. 



3. Profile Analysis 

The interstellar absorption observed toward ^ Ori is complex. Very high resolution (FWHM 
~ 0.3-1.3 km s~^) optical spectra of Na I, Ca II, Ca I, and K I reveal at least 20 components at 
heliocentric velocities between about -8 and +37 km s"^ (Welty et al. 1994, 1996, 2002; Welty & 
Hobbs 2001). These components can be associated with four low- velocity (LV) component groups, 
centered at about —2, -|-12, -|-24, and -|-37 km s^^, that are discernible in the lower resolution GHRS 
echelle and IMAPS spectra of many neutral and singly ionized species. Several weak intermediate- 
velocity (IV) components (defined here as having —60 km s"^ ^ Vq < —10 km s"^ or 45 km s~^ < 
Vq < 95 km s~^ — i.e., 27.5 km < I^lsrI ^ 77.5 km s~^), between —32 and —14 km s~^ and 
also near -|-50 km s^^, can be seen in the profiles of the stronger lines of some first ions (e.g., Mg II, 
Al II, Si II). Still weaker IV components, between —60 and —40 km s~^, may be discerned in the 
very strong C II A1334 and Si III A1206 lines. The high-velocity (HV) gas comprising Orion's Cloak 
appears at —105 km s"^ ^ v < —65 km s"^ in the strongest lines of various singly and doubly 
ionized species; the echelle profiles of those ions suggest that several components are present. The 
strongest absorption from most doubly ionized species is found in components near —16 and -|-13 
km s~^, though C III, Si III, and (probably) N III exhibit absorption over essentially the full range 
from —100 to -f-45 km s~^. Weak, broad absorption due to the higher ions C IV and Si IV extends 
from about —60 to -|-40 km s~^. Several components, which may be associated with those seen for 
the lower ions, appear to be present. A schematic diagram of the observed absorption is given in 
Figure 2. 

We have used the method of profile fitting in an attempt to discern and characterize the 
individual components giving rise to the observed absorption profiles (e.g., Welty et al. 1991, 



^Maintained at http:/ /archive. stsci.odu/copernicus/indcx. html. Sec Rogorson ot al. (1973) and subsequent articles 
for descriptions of the Copernicus (OAO-3) satellite and spectrometer and its early performance. 
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1999). In this analysis, we attempt to obtain column densities (iV), Doppler line width parameters 
(6 = [{2kT/m) + 2v^Y/'^) and heliocentric velocities {y) for as many components as are needed 
(at minimum and consistent with the measured S/N) to reproduce the observed line profiles. The 
line widths provide direct limits on the temperature (T) and turbulent velocity {yt) characterizing 
the individual components. The profile analysis provides a means for estimating the properties of 
heavily blended lines and enables us to model lower resolution spectra via the detailed structure 
derived from higher resolution line profiles (e.g., to obtain more accurate total column densities). 
We assume that each component may be represented by a (symmetric) Voigt profile (though see the 
discussion of shock models below), smeared by the appropriate instrumental response function. The 
instrumental function is assumed to be Gaussian (GHRS and IMAPS) or triangular {Copernicus), 
with the FWHM given in Table 1 or in the discussions above. We use the rest wavelengths, oscillator 
strengths, and damping constants tabulated by Morton (1991), except where more recent work has 
provided more accurate values (Table 2; see the Appendix in Welty ct al. 1999). The total sightline 
column densities, summed over all the components found in fitting the line profiles, are consistent 
with the corresponding values derived via apparent optical depth methods (e.g.. Savage &: Sembach 
1991; Jenkins 1996). 

The component structure for the high- and inter mediate- velocity gas, derived primarily from 
the high-resolution UV profiles of C II, N II, Mg 11, Si 11, and Si 111, is given in Table 3. The 
components are numbered sequentially, with the letters H and I appended to associate them with 
the high and intermediate velocity component groups, and are marked above the spectra in Figure 1. 
For the HV gas, the component parameters are most reliable for component 2H at —93 km s~^ 
(the strongest and best defined) and for component 4H at —78 km s^^ (with b constrained by the 
right-hand wing). Two additional components (IH at —97 km and 3H at —86 km s""*^) are 
required to fit the absorption shortward of component 2H and between components 2H and 4H, 
respectively; another very weak component (5H at about —68 km s~^) is found for Si 111 and C II. 
Slight asymmetries in the strongest IV absorption near —16 km s~^ (e.g., in the lines of Fe II, S II, 
and Si 11) suggest that at least two components (121 at —18 km s~^, 131 at —14 km s~^) are present 
there. IV components 61 through 91 (between —59 and —40 km s^^) are visible in the Si 111 A1206, 
N 111 A989, and C II A1334 lines. The 5-values in parentheses were determined in trial fits to the 
line profiles, but then were fixed to facilitate convergence in the final iterative fits. The component 
velocities and 6-values found in fitting the stronger HV and IV absorption features were then used 
in fitting the weaJier absorption from Al 11, Al III, S 11, S 111, Fe II, and Fe III, and also the lines of 
C III and N III observed at lower resolution with Copernicus. At least four components appear to 
be present for the higher ions C IV and Si IV, which were fitted independently. 

The absolute velocities for the various spectra were determined by fitting the low-velocity 
absorption due to (primarily) neutral gas, using the detailed component structure derived from 
very high resolution optical spectra to model the lower resolution UV spectra (Welty ct al. 1999). 
The measured velocity offsets between the GHRS echelle spectra and the optical spectra ranged 
from about +A km s~^ for the first ECH-A spectra obtained to about —3.5 km s~^ for the last 
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ECH-B spectra. The corresponding offsets for the IMAPS and Copernicus spectra are about 2 km 
s~^ and 7 km s~^, respectively. 

4. Abundances, Depletions, and Physical Conditions 

Table 4 gives total column densities (or limits) for various neutral and ionized species in the 
HV and IV gas. The values for the heavy element species are based on the profile fits and/or the 
equivalent widths; the estimates for A^(H I) and A/'(H II) are discussed below. Figure 3 shows the 
HV and IV column densities as functions of ionization potential. By combining near-UV spectra 
from GHRS with far-UV spectra from IMAPS and Copernicus, we have HV and IV column densities 
for both singly and doubly ionized states of C, N, Al, Si, S, and Fc. We also have significantly 
restrictive limits for the corresponding neutral species and (for the HV gas) for the triply ionized 
states of C, Si, and S; N IV is also unlikely to be significant in the HV gas, given the limits on C IV, 
Si IV and S IV. 

The elements nitrogen and sulfur typically have nearly solar abundances (within about 20%) 
in diff'use Galactic clouds; carbon is generally very slightly depleted (by a factor of about 2.5); 
Al, Si, and Fe are usually more severely depleted (factors of 4-300) (e.g., Savage &; Sembach 
1996; Fitzpatrick 1996; Welty et al. 1999; and references therein). We use the meteoritic (CI 
carbonaceous chondritic) reference abundances of Anders & Grcvcsse (1989) for all elements except 
C, N, and O, for which we use solar photosphcric values from Grcvesse &: Noels (1993).^ The total 
hydrogen column densities A^(H l)+A^(H II) for the HV and IV gas may therefore be estimated 
from the sums iV(N II)+7V(N III), A^(S Il)+iV(S III), and (perhaps) iV(C Il)+iV(C III); the level of 
depletion may be gauged via the corresponding sums for Al, Si, and Fe (assuming negligible Al IV 
and Fe IV) . The column density limits for O I, the dominant ionization state of oxygen in neutral 
(H I) gas, may be used to set corresponding limits on the column density of H I associated with 
each component group. 

Table 5 summarizes some of the physical properties of the individual HV and IV components 
toward C Ori. Using estimates for the temperature derived from the component line widths (§ 4.1), 
analysis of the C II fine-structure excitation equilibrium yields estimates for the electron density 
(ne), which is approximately equal to the local hydrogen density (nn) for these predominantly 
ionized components. (If helium, with first ionization potential 24.6 eV, were entirely singly ionized, 
then He would be roughly 10% higher than nn-) Dividing the total hydrogen column densities by 
nn then yields the thicknesses (r) of the individual components. Table 6 compares various ion 
ratios and derived physical properties (averaged over all components) for the HV and IV gas seen 
toward C, Ori, 23 Ori (Welty et al. 1999), and r CMa (appendix) with predictions from shock 



^Use of the revised solar photospheric abundances for the few elements reported by Holweger (2001) would have 
a minor impact on the results discussed in this paper. 
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models (discussed below). 

4.1. High- Velocity Gas (—105 km ^ ^ —65 km s~^) 

The column densities of C II + C III, N II + N III, and S II + S III suggest that C, N, and S are 

present in nearly solar proportions (within ±25%) in the HV gas — and thus that carbon (as well as 
nitrogen and sulfur) may be at most only very mildly depleted. Under the assumption that nitrogen 
and sulfur are undepleted, we estimate log[A^(H I) + A^(H II)] (cm~^) ~ 17.88±0.06 for the HV gas. 
Consideration of the uncertainties on A^(S II) and A^(S III) and the upper limit on N{S IV) yields 
the fairly strict upper limit log[iV(H I) + iV(H II)] (cm~^) < 17.96. The uncertainties in the total 
iV(H) (and in the resulting depletions of C, N, and S) are dominated by the uncertainties in A^(C III) 
and A"(N III) (both derived from strong lines in lower resolution spectra) and in the adopted solar 
reference abundances for the three elements (0.05-0.07 dex). Carbon may be very slightly depleted 
(by roughly 0.1 dex) in the HV gas, but within the uncertainties could be undepleted (and we note 
that the ratio of carbon to sulfur is essentially solar). 

We also have fairly stringent upper limits on high-velocity absorption from N I and (especially) 
O I; the limit log[iV(H I)] (cm^^) < 14.84 in Table 4 assumes that oxygen is undepleted. The HV 
gas toward ( Ori is thus primarily ionized, with iV(H I)/7V(H II) < 0.001. The lack of high- velocity 
absorption from the higher ions S IV, Si IV and C IV suggests that only a narrow range of ionization 
states is significantly populated for any element. 

The HV column densities for Al II + Al III, Si II + Si III, and Fe II -|- Fe III all correspond to 
somewhat smaller H II column densities, suggesting some depletion of Al (—0.82 dex), Si (—0.39 
dex), and Fe (—0.97 dex), relative to C, N, and S (see discussion in § 5.4). Since C III and N III 
are present (and dominant) in the HV gas, there could be some Al IV and Fe IV as well, however, 
so that the depletions of aluminum and iron could be less severe. For temperatures consistent with 
the observed ion ratios (see below), coUisional ionization equilibrium models (e.g., Sutherland & 
Dopita 1993) suggest that Al IV may well be significant, but not Fe IV. The shock models described 
below, however, suggest that Fe IV might also be significant. 

Given the definition for the line width parameter b, we may obtain estimates for both the 
temperature (T) and internal turbulent velocity {vt) of the individual HV components by comparing 
the 6- values for absorption lines due to elements of different atomic weight. For component 2H, 
where the line widths are fairly well-determined for C II, N II, Mg II, Si II, and Fe II, the individual 
6- values listed in Table 3 yield firm upper limits Tmax ^ 12,000 K in most cases, and the ensemble of 
6-values is most consistent with T ~ 9000ib2000 K and vt ~ 0.5ib0.5 km s~^. Since the 6-values for 
the other HV components are not as well determined, we will assume that all the HV components 
have T ~ 9000 K. Such a temperature is significantly smaller than would be inferred from the 
various ion ratios listed in Table 6, if the HV gas were assumed to be in CIE. The HV gas is 
therefore likely not in CIE, but appears to have cooled from an initially higher temperature more 
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rapidly than it could recombine, perhaps after having been shocked; photoionization may also be 
significant (see discussions below). As the values for the HV Si III components can be larger than 
those for the corresponding Si II components, at least some fraction of the higher ions may be at 
somewhat higher T and/or more broadly distributed. 

In this ionized HV gas, the excited state C II* is populated primarily via collisions with elec- 
trons, and the ratio A''(C II*)/A''(C II) can be used to estimate (and thus nn) from the equation 
of excitation equilibrium, since the temperature is known. In the limit of small A'"(C II*)/iV(C II), 
rie = 0.186 [N{C I1*)/N{C II)] (Trapero et al. 1996; Welty et al. 1999). For T ~ 9000 K, 
the He ^ nn derived for HV components 2H and 4H are both ~ 0.1-0.2 cm~^. The limits for 
components IH and 5H are consistent with such values, but the upper limit for component 3H is 
a factor 2-4 smaller (Table 5). The thermal pressures for components 2H and 4H, as measured 
by 2ncT, are thus about 2000-4000 cm~^K. Estimating individual HV component A'^(H II) from 
A'"(S II) + A^(S III), and adjusting for the adopted shghtly higher total A'"(H II), the thicknesses of 
components 2H and 4H are each about 0.5-0.6 pc; component 3H appears to be somewhat thicker. 

4.2. Intermediate- Velocity Gas (—60 km ^ vq < —10 km s~^) 

For a number of the singly and doubly ionized species, the column densities in the IV gas are 
typically about a factor of 3 higher than in the HV gas (Table 4). From N{S II) + iV(S III) [and 
the limits on N{S I) and N{S IV)], we estimate that log[iV(H I) + iV(H II)] (cm'^) ~ 18.44 for 
the IV components. The corresponding sums for carbon and nitrogen are not as well determined 
(the C III A977 line is very strong; the N III A989 line is strong and blended with a Si II line), but 
suggest very similar IV N{IV) (and negligible depletion for carbon and nitrogen). 

In contrast to the HV gas, some O I is detected in the strongest IV components near —16 
km s~^. The implied column density log[A^(H I)] (cm^^) ~ 15.31 then yields a neutral fraction 
iV(H I)/iV(H II) ~ 0.001 for the IV gas, so that the IV gas is also predominantly ionized. The 
N{X III)/iV(X II) ratios are generally within a factor of 2 of the corresponding ratios in the HV 
gas (Table 6), suggesting that the degree of ionization is similar for both. 

The IV column densities for Al II + Al III suggest that aluminum may be depleted by about 
—0.86 dex, relative to sulfur, in the IV gas close to the aluminum depletion inferred for the HV 
gas. Unfortunately, the IV column densities for Si III and Fe III are not well determined — the 
Si III A1206 line is saturated over part of the IV range and the Fe III A1122 line profile is noisy and 
incomplete — so that we do not have strong constraints on the depletions of silicon and iron in the 
IV gas. As was noted for the HV gas, the depletion of aluminum in the IV gas would be less severe 
if some Al IV is present. 

The 6-values adopted for most of the IV C II, Mg II, Al II, Si II, and Fe II components suggest 
that the temperatures are comparable to those in the HV gas (i.e., about 9000 K). Analysis of 
the C II fine-structure excitation, assuming T ~ 9000 K, yields He ~ 0.16 cm~^, and thus 2neT ~ 
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2500-3000 cm~^K, for several of the stronger IV components — very similar to the densities and 
thermal pressures found for the HV gas (Table 5). The IV components have total A^(H II) higher 
by a factor of about 3, however, and so the total thickness (~ 5.1 pc) is correspondingly larger. 

While the lower ionization stages are generally dominant in the IV gas, some of the absorp- 
tion from the higher ions S IV, Si IV, C IV, and (perhaps) N V may be associated. The ratio 
A^(C IV)/A'^(Si IV) = 7.0ib0.9 found for the IV gas is significantly higher than the median value 
3.8±1.9 for 31 lines of sight observed with lUE and/or GHRS (Sembach, Savage, k Tripp 1997), 
but is similar to the ratio found at intermediate velocities toward 23 Ori (Welty et al. 1999). The 
GHRS echelle spectra suggest that the broad, weak absorption from C IV and Si IV toward ( Ori 
arises from several narrower components, some at velocities similar to those found for components 
in lines of lower ionization stages such as Al III, S III, and Si III (e.g.. Savage, Sembach, & Cardelli 
1994; Fitzpatrick & Spitzer 1997). 



5. Discussion 

5.1. Orion's Cloak and Other High- Velocity Gas in the Galactic Disk 

Comparisons among the 21 cm. Ha, soft X-ray, and 100 //m emission maps of the Orion- 
Eridanus region have suggested the presence of a large, expanding shell, predominantly neutral but 
with an inner ionized zone, encompassing a bubble of hot gas at T ~ lO^-^ K (Heiles 1976; Reynolds 
& Ogden 1979; Burrows ct al 1993; Brown et al. 1995; Snowden et al. 1995). Cowie et al. (1979) 
ascribed the high-velocity gas seen with Copernicus toward a number of Orion stars (including ( 
Ori) to a radiative shock due to one or more recent supernovae in that region. They estimated a 
radius of about 120 pc (i.e., outside the expanding neutral shell) and an age of about 3 x 10^ yr for 
this apparently large scale feature, as well as a pre-shock ambient density of about 3 x 10^'^ cm^^. 
Ionized gas at intermediate velocities is also seen toward a number of the Orion stars, suggesting 
that the IV gas also is a large-scale, regional feature — not intimately associated with the individual 
stars. The HV and IV gas is not seen toward all the stars in the region, however, suggesting that 
its distribution is somewhat patchy (see also Shore 1982; Sonneborn et al. 1988). Cowie et al. 
found the IV gas to be more highly ionized than could be attributed to collisional ionization at 
those velocities, and concluded that the IV gas is photoionized by (but at some distance from) the 
Orion association stars. 

Higher resolution GHRS and IMAPS spectra covering additional species have enabled a more 
detailed understanding of the composition and ionization of this high- and inter mediate- velocity 
gas, and have yielded estimates for the temperatures and electron densities in the individual HV 
and IV components (Trapero et al. 1996; Welty et al. 1999; this paper). In Table 6, we compare 
total hydrogen column densities, ion ratios, depletions, and average physical properties for the 
HV and IV gas seen toward ^ Ori, 23 Ori, and r CMa (the latter somewhat removed from the 
Orion's Cloak region; see the Appendix). In all three sightlines, the high- velocity gas appears to 
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consist of 4-5 components, whose similar widths (for a given ion) suggest that the HV component 
temperatures are typically in the range 7000-11000 K in all three cases. The total HV H II column 
densities are comparable toward 23 Ori and ( Ori, but A^(H II) is smaller by a factor of about 
6-7 toward r CMa. The available N{X III)/A^(X II) ratios suggest that the HV gas toward C Ori 
(09.5 lb) is more highly ionized than the HV gas toward 23 Ori (Bl V), which in turn is more 
ionized than the HV gas toward r CMa (09 lb) — so that the degree of ionization is not correlated 
with the spectral type of the background star. The depletions of silicon and aluminum in the HV 
gas (relative to typically undepleted sulfur) arc slightly more severe toward ( Ori than toward 23 
Ori, and arc least severe toward r CMa. Carbon, which is typically depleted by about 0.4 dex in 
both diffuse and moderately dense disk clouds at lower velocities (Sofia et al. 1997), appears to be 
essentially undepleted in the HV gas toward all three stars. The density and thermal pressure are 
higher by a factor of about 3 in the HV gas toward 23 Ori, compared to the other two sightlines. 

The total intermediate-velocity H II column density is about 3 times higher toward ( Ori than 
toward 23 Ori; no IV components were detected toward r CMa [except, perhaps, in Copernicus 
spectra of C III (Cohn & York 1977)]. Comparisons of the A^(X III)/Ar(X II) ratios between the 
HV and IV gas, toward both ( Ori and 23 Ori, suggest that the ratios may be slightly higher in 
the IV gas (though the ratio for carbon may be slightly lower in the IV gas toward ( Ori). The 
depletions in the IV gas are similar to those found for the respective HV components toward both 
C Ori and 23 Ori (with less severe depletions toward 23 Ori). The IV components toward 23 Ori 
have significantly higher densities and thermal pressures than those toward C, Ori. 

We may also compare the high- and intermediate-velocity gas observed toward ( Ori and 23 
Ori with gas found at comparable velocites in other regions. For example, absorption at |vlsr| ^ 60 
km s~^ (both positive and negative velocities) has been observed toward a number of stars located 
behind the Vela supernova remnant. Jenkins, Silk, & Wallerstein (1976) analyzed Copernicus 
spectra of four stars, and found C, N, Si, and Fe to be essentially undepleted in the high-velocity 
gas at —180 km s^^ toward HD 74455 and at —90 km s~^ toward HD 75821. Limited data for 
various ratios N{X IIl)/iV(X II) suggest that the HV gas toward Vela is not as ionized as the 
HV gas toward Orion. Jenkins, Wallerstein, & Silk (1984) examined lUE spectra of 45 stars in 
the Vela region, and found the depletions of Al, Si, and Fe in the HV gas to be less severe than 
— 1.2 dex, —0.5 dex, and —0.55 dex, respectively (relative to S) — and thus less severe than the 
corresponding depletions in lower velocity clouds in the same region ( — 2.2 dex, —0.65 dex, and —1.2 
dex). Raymond, Wallerstein, &: Balick (1991) compared UV absorption lines toward HD 72088 with 
emission lines from adjacent nebulosity, attributed the HV gas to a shock with Vg ~ 150 km s~^ 
and a swept up column density of about 10^^ cm~^ (comparable to the values found for the HV 
gas toward Orion), and surmised that carbon might be undepleted. Jenkins et al. (1998), using 
HST/STIS spectra, compared the column densities of various species normally dominant in H I 
regions for the gas at —70 km s^"*^ < Va^ < —45 km and 112 km s^^ ^ o ^ 130 km s^^ toward 
HD 72089. The low A^(0 l)/iV(S II) and iV(N I)/iV(S II) ratios suggest that the gas is ionized 
(though the ratios are not quite as low as those in the HV gas toward Orion). The depletions 
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of Al, Si, Fe, and Ni appear to be very mild (—0.6 dex, —0.3 dex, —0.4 to —0.8 dex, and —0.5 

dex, relative to sulfur) — comparable to the values found in the HV gas toward Orion [though 
A^(Fe II)/A^(S II) is lower toward Orion, where Fe III appears to be the dominant ionization state 
of iron]. Determination of various ion ratios and physical properties characterizing the HV gas 
toward the Vela SNR would enable more detailed comparisons with the HV gas toward Orion; new 
data are required. 

5.2. Ionization Corrections / Models of Ionized Gas 

Determinations of interstellar abundances and depletions are often based on observations of 
a single ionization state — usually presumed to be the dominant one for the gas in question. In 
predominantly neutral (H I) gas, the first ion (X II) is dominant for most of the more abundant 
elements (those with first ionization potential less than about 13.6 eV), so that depletions may be 
fairly accurately estimated via ratios such as A^(X II)/A^(S II) or N(K H)/N(Zn II), since both 
sulfur and zinc are typically nearly undepleted (e.g., Welty et al. 1997). Such simple estimates may 
not yield accurate relative abundances or depletions, however, for lower column density clouds (e.g., 
A''(H) < 10^^ cm~^), for significantly ionized gas, and for clouds where the degree of ionization is 
not known (e.g., Sembach et al. 2000; Howk & Sembach 1999). Unfortunately, data for potentially 
significant higher ions — especially at a spectral resolution high enough to compare component 
structures with those found for the lower ions — often are not available. Of the elements we have 
considered in this paper, for example, only aluminum has detectable lines longward of Lyman a 
from its second ion (Al III). While lines from Si III and S III just shortward of Lyman a can be 
observed with HST, the permitted resonance lines for C III, N II, N III, and S IV all lie below 1100 
A. For studies of the Galactic ISM, IMAPS spectra are available for a small number of fairly bright 
stars; otherwise, we must depend on lower resolution spectra obtained with Copernicus or with 
FUSE. For studies of QSO absorption-line systems, the shorter wavelength lines from those higher 
ions are entangled in the Lyman a forest. Ionization corrections therefore are often estimated from 
theoretical models, by assuming a particular cloud structure (including coexistence of various ions) , 
a radiation field (or lack thereof), and either coUisional or photoionization equilibrium. 

Because we have column densities for both singly and doubly ionized C, N, Al, Si, S, and 
Fe toward ^ Ori, we may (1) estimate empirical ionization corrections for the HV gas and (2) 
make more extensive comparisons with the predictions of various models of ionized gas. For all six 
elements, the doubly ionized form is dominant in the HV gas; the ratio A^(X III)/A^(X II) ranges 
from about 1.4 for nitrogen to about 11 for iron (Table 6). In Table 7, we compare the observed HV 
gas ion ratios N{X ll)/N{S II) and A''(X m)/N{S III) with the ratios of the total column densities 
N{X)/N{S) — all relative to the corresponding solar ratios listed in column 2. Because the second 
ions are dominant in the HV gas, the [X III/S III] ratios provide the best estimates for the overall 
abundance ratios; all are within ±0.1 dex of the corresponding [X/S].^ In this case, however, the 



[X II/S II] ratios do nearly as well — within ±0.2 dex — except for iron. 



In Figure 4, we compare the observed HV gas ratios iV(X III)/Ar(X II) and iV(X IV)/iV(X II) 
(given by the filled circles in each panel) with the values predicted by several models for ionized 
gas: 

• Sembach et al. (2000) presented ion fractions for models of relatively low density photoionized 
gas calculated using the CLOUDY code (e.g., Ferland et al. 1998). In those models, Ti, Cr, 
Mn, and Fe can be doubly ionized more readily than most other elements, so that the ratios 
N{X II)/iV(S II) yield underestimates of the relative total abundances X/S for those four 
elements. For most other elements, the ratios N{X lI)/N{Fe II) correspondingly overestimate 
the true X/Fe. The magnitudes of the differences depend somewhat on the shape of the 
radiation field, but can easily exceed factors of 2. The HV gas toward ( Ori is significantly 
more ionized than the gas in the Sembach ct al. models, however, as the observed ratios 
N{X lll)/N{X II) (for C, N, Al, Si, S, and Fe) are all more than an order of magnitude larger 
than the values predicted for their "standard composite" model (open circles in top panel). 

• Gruenwald &: Viegas (1992) tabulated ion column densities for lines of sight through model 
H II regions characterized by several representative densities, metallicities, and radiation 
fields. For the models with solar metallicity, nn = 10 cm~^ (the lowest density consid- 
ered), and an 04 star radiation field (open squares in top panel), most elements have 
N{X lll)/N{X II) > 1.5 near the center (upper point), decreasing to 1 by a fractional 
radius r of 0.7 (lower point); sulfur is systematically more ionized, however. While the ion 
ratios for carbon and nitrogen are (at some radii) similar to those seen in the HV gas toward 
( Ori, the ratios for silicon and iron are lower and the ratio for sulfur is higher. For r = 
0.7, the predicted N{X TV)/N{X II) are consistent with the limits found for the HV gas — 
but are much larger than those limits for smaller r. We note, however, that 0^ Ori A and 
C appear to have the earliest spectral types (07 and 06) of any of the stars in the Orion 
region; most of the brighter stars in the region have types between 09 and Bl (e.g., Warren 
(fe Hesser 1977; Wall et al. 1996). It is therefore unlikely that the radiation field toward C Ori 
could be characterized by as early a spectral type as 04. For the softer BO star field (open 
triangles in top panel; single points are for r = 0.0), the first ions are generally dominant 
throughout (again except for sulfur); the ion ratios are thus lower for all r and all elements 
(except sulfur) than those seen in the HV gas. The ratios predicted for a slightly harder 
field (e.g., 09) and r ^ 0.0 might be reasonably consistent with most of the observed ratios 
— except for iV(Fe IIl)/iV(Fe II) (too low) and N{Si IV)/iV(Si II) (too high). Comparisons 
among the various models suggest that the N{X III)/iV(X II) ratios would be somewhat lower 
for densities below 10 cm~^. 



The square braces denote abundance ratios relative to the adopted solar ratios: [X/Y] = log(X/Y) — log(X/Y)0. 
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• Sutherland Sz Dopita (1993) have computed models for gas in collisional ionization equilib- 
rium. The middle panel of Figure 4 compares the ion ratios predicted for CIE at T = 20,000 
K (open squares), 40,000 K (open triangles), and 80,000 K (open circles) with those in the 
HV gas (filled circles). The CIE models are consistent with the observed ionization of C, N, 
and O for T ~ 50,000-80,000 K, and with the observed ionization of Mg, Al, Si, S, and Fe 
for T ~ 25,000-40,000 K. As noted above, however, the observed line widths imply a much 
lower temperature (near 9000 K), so the HV gas must not be in CIE. 

• Even models of cooling gas (e.g., Shapiro Sz Moore 1976; Schmutzler Sz Tscharnuter 1993; bot- 
tom panel, open squares and open triangles, respectively), however, predict lower N{X 111)/N{X II) 
ratios at T = 10^ K than those observed in the HV gas — perhaps reflecting the neglect of 
ionization via radiation (stellar and/or shock-related) in the models. On the other hand, the 
ratios N{C IV)/iV(C II) and N{Si IV)/iV(Si II) at T = 10^ K predicted by Shapiro & Moore 
(1976) are much higher than the observed upper limits — though comparison with the results 

of Schmutzler Sz Tscharnuter suggests that those predicted ratios will decline rapidly below 
T = 10^ K. 

Comparisons between the observed HV gas ion ratios and the predictions from models of 
diffuse photoionized gas, denser H II regions, collisionally ionized gas (in equilibrium), isochorically 
cooling gas, and shocked gas (see next section) thus indicate that a wide range of ion ratios may be 
produced in those different physical situations — dependent both on current physical conditions 
and on the history of the gas. If sufficient diagnostics are available, those differences can lead 
either to an accurate understanding of the observed gas or to the realization that other (as yet 
unexplored) processes must be significant. Of the models noted above, photoionization by an 09 
star radiation field might yield the best match with the observed ion ratios [based on a crude 
interpolation between the results for the 04 and BO star models of Grucnwald & Vicgas (1992)] 
— if the field due to the ensemble of stars in the region can mimic that computed for r ~ 0.0, 
and if the much lower observed density (0.1 cm~^ versus 10 cm~^) does not significantly affect the 
predicted ratios. While some of the models predict some ion ratios consistent with those observed 
for the HV gas, none of the models appears able to reproduce all of the observed ion ratios. 

5.3. Shock Models for High- and Intermediate- Velocity Gas 

The presence of ionized gas at high and intermediate velocities relative to the LSR; the prox- 
imity of that gas to the numerous young, massive stars in Orion; evidence from other data for an 
expanding shell produced by supernovae and/or stellar winds; the relatively mild depletions of Si, 
Al, and Fe; and indications that the gas is not in CIE all suggest that shocks may have played 
a major role in the evolution of that gas. The high-resolution, high S/N UV spectra of ( Ori 
have yielded information on the line profiles, ion ratios, and physical conditions in the HV and 
IV gas which can provide good tests for theoretical models of interstellar shocks. In particular. 
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comparisons of predicted and observed ratios of different ionization states of a single element [e.g., 
A^(Si IIl)/iV(Si II)] avoid possible complications due to differences in depletions (interesting in their 
own right) which may arise in comparisons of ions of different elements [e.g., A^(C IV)/A^(Si IV)]. 
We have considered two sets of shock models to try to understand the origins of the HV and IV gas 
toward ( Ori: (1) for the HV gas, a shock moving toward us at a velocity of order —100 km — 
corresponding to the Cowie et al. (1979) model, and (2) for the IV gas, a shock moving away from 
us, produced when a flow of high-velocity gas encounters lower velocity material — as proposed 
by Jenkins & Pcimbcrt (1996) to explain the systematic broadening and small velocity shifts, with 
increasing rotational level J, of the H2 lines near w ~ km s~^ seen in the IMAPS spectra of ^ 
Ori. In that second model, the HV gas is the pre-shock flow. 



5.3.1. Models for the High- Velocity Gas 

Ion column densities and absorption line profiles were computed for radiative shock wave 
models using an updated version of the computer code described by Raymond (1979) and Cox & 
Raymond (1985). In brief, the model parameters are shock speed (vg); pre-shock density (no), ion- 
ization state, and magnetic field (Bq); and elemental abundances. Because of numerical limitations 
in the code, the pre-shock densities are taken to be an order of magnitude or more larger than 

would produce ^ rin = 0.1-0.2 cm^"' in the C II region (as derived above). As long as we scale 

1/2 

the magnetic field and the radiation field appropriately (as and no, respectively), however, 
this will have no effect at all on the resulting column densities or line profiles. The shock jump 
conditions (e.g. McKee &; Hollenbach 1980) are used to find the density, temperature and magnetic 
field strength just behind the shock. The code then follows the time-dependent ionization and 
radiative cooling of the shock-heated gas; photoionization by radiation emitted during the cooling 
is included. In some of the calculations presented here, we also include ionizing radiation (VFrad) 
from the Orion stars, approximated by a diluted Kurucz (1979) 25,000 K model spectrum. All the 
models used the abundances of Baldwin et al. (1991). 

Models for absorption lines require an additional model parameter that can often be ignored 
in emission line models. When the shocked gas cools to about 1000 K, it ceases to emit signifi- 
cant radiation at optical and UV wavelengths, and it does not matter how much cooler material 

accumulates. For absorption line models, however, the integration must be cut off at some point. 
Fortunately, the cutoff column density, log N(H) (cm~^) = 17.88, is tightly constrained for the C 
Ori high-velocity component. 

For a strong shock, the flow speed is 3/4 Vg just behind the shock, and approaches Vg as the 
cooling gas is compressed. At each point in the flow, we compute the opacities in some of the 
strong UV lines, assuming a thermal velocity distribution at the local temperature, shifted to the 
local bulk velocity. These contributions to the optical depth are summed to produce the predicted 
absorption line profiles for the shock along a line normal to the shock front. Viewing the flow at an 
angle from the shock normal will increase the column densities by a factor of sec 6 and compress 
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the velocity scale by a factor of cos 9. 

In Figure 5, the observed HV profiles are compared with those predicted by several models, 
whose parameters are listed in Table 8. Some predicted ion ratios are given in column 5 of Table 6 
and arc plotted in Figure 6. Model A (dotted lines) represents a "simple" 95 km shock; models 
B (short-dashed lines) and C (long-dashed lines) add (in succession) a stronger radiation field and 
a stronger magnetic field. Increasing the stellar radiation field (model B) generally decreases the 
column densities of the various singly ionized species and enhances the higher ions. Increasing the 
magnetic field (model C) has somewhat the opposite effect — the column densities of the singly 
ionized species are enhanced, while those of C IV and (to a lesser extent) Si IV are reduced. A 
magnetic field also limits the compression of the shocked gas and reduces the velocity offset between 
high- and low-ionization species. 

The observed HV gas line profiles appear to be in qualitative agreement with the basic shock 
picture. A 95 km s~^ shock coming towards the Earth will produce narrow absorption features 
at that speed in the lowest ions (see, e.g., the observed C II, Si II, and Mg II profiles), with an 
asymmetric extension to higher (i.e., less negative) velocities — as well as C III, N III, Si III, and 
S III column densities similar to those observed. The much weaker lines from higher ionization 
stages (e.g., C IV, Si IV) are predicted to be less blue-shifted by up to Vg/^ ~ 24 km s~^. Because 
the ionization lags behind the temperature and velocity evolution of the shock, the shifts for the 
higher ions in these models are less than half that maximum value. 

There are several major differences between the observed and predicted profiles and ion ratios, 
however: 

• First, the predicted blue shift is not quite large enough to match the observations unless the 
shock speed exceeds 100 km s~^ — but such fast shocks produce too much C IV and Si IV 
absorption. That problem is easily solved, as we have no reason to believe that the shock 
is moving through stationary gas, and random velocities of 10-20 km s~^ are common in 
the ISM. We therefore choose the speed for the pre-shock medium {vq) that shifts the strong 
absorption components to ~ —95 km s~^. Alternatively, adding a strong pre-shock magnetic 
field would absorb some of the shock energy — thus permitting a higher speed shock without 
creating higher ionization. A strong field would also reduce the shock compression, however, 
so that the predicted speeds would not increase. 

• Second, while model A predicts N{X. III)/A?^(X II) ratios for C, N, and S that are lower than 
those observed (by factors of order 2) , the corresponding predicted ratios for silicon and iron 
are too high by factors of 4-5 (Figure 6). The disagreements between observed and predicted 
ion ratios are typically somewhat worse for models B and C. All three models predict too 
much Si IV; models A and B predict too much C IV. It is perhaps interesting that the models 
predict too much highly ionized material for the refractory elements silicon and iron. Such 
an effect might be expected if a significant fraction of the observed material is liberated from 
grains by the shock. The refractory elements are only gradually (but increasingly) sputtered 
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from the grains as the gas cools (e.g., Jones et al. 1994) — so they are still severely depleted in 
the highest temperature gas. Therefore, (1) only a small fraction of those atoms would reach 
the highest ionization states predicted by the models, and (2) the atoms will be returned to 
the gas from the dust predominantly as lower ions. For elements that are initially severely 
depleted, the return of even 10% of the atoms to the gas phase (see § 5.4 below) could have 
noticeable effects on both the ion ratios and the line profiles. 

• Third, the profile widths in the model spectra are systematically of order 50% broader than 
those of the strongest observed absorption component 2H (note that the model profiles have 
not been convolved with the relevant instrumental response functions in the figure). This 

discrepancy may suggest that the post-shock gas has not cooled enough in the models by the 
time the total column density has reached the cut-off value. While the difference in widths 
also could be due to observing the shock at an angle of about 45 degrees, the large scale 
nature of the HV gas in Orion would seem to suggest that the line of sight to Ori should be 
roughly normal to the shock. Furthermore, the corresponding increase in shock speed would 
produce too much C IV and Si IV. 

• Fourth, the models produce a tail toward lower speeds in the profiles of the singly ionized 
species, but it is far narrower than the 20 km s~^ tails in the observed C II and N II profiles. 
In the HV gas toward 23 Ori (Welty et al. 1999b) and r CMa (appendix), the corresponding 

redward absorption appears to be due to distinct additional components, as we have assumed 
in fitting the HV gas profiles toward C Ori. (We note, however, that the models indicate 
that shocks can produce noticeably asymmetric profiles, contrary to what we have assumed 
in fitting the observed spectra.) 

A number of interesting physical effects might contribute to the extended lower speed tails. 
Inhomogeneities in the density or magnetic field on a length scale smaller than the cooling length 

(which is enormous, given the low density) are probably most likely. Such inhomogeneities could 
affect the shock speed (rig) or the shock compression (Bq). Detailed computation of these effects is 
beyond the capability of the 1-D steady flow models presented here, but we estimate that density 
fluctuations of 60% could account for the observed overall 30 km spread in component speeds. 
The required magnetic field fluctuations depend on the values of no and Bq, but for the no in Model 
A, small regions with 5o ^ 1 /xG could provide the lower speed gas. 

The low densities derived from the N{C ll*)/N{C II) ratios in the HV components — rie ~ 
nn ^ 0.1-0.2 cm~^ — suggest that (1) the pre-shock densities must have been very low (perhaps 
of order 0.001 cm~^), and (2) the shock may be better characterized as a slow mode (rather than 
a fast mode) MHD shock (e.g., Draine & McKee 1993). The radiative cooling time just behind 
the shock is of order 1500/no years — which can be longer than the expansion time for adiabatic 
cooling for such low pre-shock densities. Again, computation of detailed line profiles for such shocks 
is beyond the capabilities of the models considered here. 
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5.3.2. Models for the Intermediate-Velocity Gas 

We have also compared the absorption-Une profiles of the intermediate- velocity gas with shock 
models. The presence of weak C IV and Si IV, along with an offset of about 10 km s~^ between 
the low and high ionization species, suggests a shock speed of about 90-100 km s~^ (similar to 
the values explored for the models of the HV gas). Faster shocks produce too much C IV, Si IV, 
and S IV, while slower shocks produce none. The various ion ratios A^(X III)/A^(X II) and the 
(relatively low) density found for the IV gas also are similar to those found for the HV gas. 

The IV gas differs from the HV, however, in that its column density is larger by about a factor 
of 3, its velocity is only a quarter that of the HV gas, and the sense of the offset between high and 
low ionization states is inverted. The latter two aspects suggest that the direction of the shock is 
reversed compared to that of the HV gas. If gas moving at 100 km s~^ strikes a stationary dense 
cloud, the gas just behind the shock is compressed by a factor q and moves toward the cloud at a 
speed (100/g) km s~^. As the gas cools and q increases, its speed decreases. 

In Figure 7, we compare the predicted profiles for the models with initial flow velocities of —95 
to —100 km s~"^ (Table 8) with those observed for the IV gas toward C, Ori; the predicted ion ratios 
are given in column 8 of Table 6 and are shown in Figure 8. Model D (dotted lines) is the "basic" 
one here; models E (short-dashed lines) and F (long-dashed lines) add (in succession) a stronger 
magnetic field and a stronger radiation field. 

As for the HV gas, there is qiialitative agreement between the predicted and observed line 
profiles. In particular, for these shocks moving away from us, the asymmetric tail is now toward 
more negative velocities, with an apparent progression in the profile shapes from high to low ions 
(e.g., for N III, N II, and N I in Fig. la). The higher ions have more prominent absorption in the 
asymmetric tail, while the lower ions are more concentrated near ~ km s~^ — as would be 
expected for post-shock gas that starts hot at the most negative velocity, then cools and collapses 
toward « ~ km s~^. Similarly, there also appears to be a systematic increase with v in ngT""^/^, 
as indicated by the ratios N{C II*)/iV(C II) and iV(N II*)/A^(N II) — consistent with cooling, 
compressing gas. 

Again, however, there are also significant differences between the observed and predicted pro- 
files: 

• First, the predicted absorption shortward of about —30 km s~^ in the strong lines of C II, 
C III, and Si III is generally weaker than the observed absorption. A lower compression in the 
shock (e.g., by a factor 2.6 instead of the canonical factor 4) — perhaps due to a magnetic 
field — might help to enhance the absorption at more negative velocities (Jenkins &; Peimbert 
1996). As for the HV gas, the apparent discrete components seen in the IV absorption for 
several species — which are not seen in the smooth predicted profiles — could be due to 
inhomogeneities in the pre-shock density and/or magnetic field. 
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• Second, the predicted ion ratios do not agree with the observed ratios (Fig. 8). For models 

D and E, the predicted column densities of most singly ionized species are somewhat higher 
than the observed values (especially for silicon and magnesium), while the corresponding 
iV(X III) are lower than the observed values (generally by factors of at least 4). Moreover, 
these two models predict fairly strong IV absorption from N I, which is not observed. The low 
N{X III)/A/"(X II) ratios and the strong N I both are likely due to having no stellar radiation 
field in those two models. For model F, which includes a strong stellar field, most of the ion 
ratios are higher than those observed for the IV gas. A somewhat weaker radiation field than 
that adopted in model F — perhaps similar to that in model A — might yield better overall 
agreement with the observed ion ratios. 

• Third, while models D and E predict total IV column densities for C IV and Si IV close to 
the observed values, the observed absorption covers a wider velocity range. Model F predicts 
somewhat higher IV column densities for those higher ions. 



The most serious problem with the reverse shock scenario, however, is the nature of the pre- 
shock gas, which should be visible near —100 km s~^ — unless (by a remarkable coincidence) the 
shock happens to have passed through most of the high-speed gas. One could imagine that the 
HV components are just clouds of gas expanding away from the Orion stars (rather than a shock), 
and that the IV gas is the post-shock flow following a reverse shock. However, a shock capable 
of slowing and heating the gas as much as is required would also compress the gas by perhaps an 
order of magnitude — but the A''(C ll*)/N{C II) ratios suggest essentially equal densities in the 
HV and IV gas. In view of this and the other problems listed above, the IV absorption may not 
result from a shock such as those described in this paper. The more complex switch-on, switch-off 
shock scenario (with compression as low as a factor 2.6) proposed by Jenkins &; Peimbert (1996) 
is, unfortunately, beyond the capabilities of these models. 

The abundances and physical properties characterizing the gas near f ~ km s~^ (the "low- 
velocity obstruction" in this reverse shock scenario) are also of interest. Fits to the proflles of Zn II 
and S II suggest that the components near f ~ km s~^ contain about 10% of the neutral gas 
toward ( Ori. Weak absorption from N II, Si III, and S III indicate that a small amount of ionized 
gas is present as well [with iV(X III)/iV(X II) < 0.01 for X = Al, Si, and S]. Column densities for 
Si II, Al II, and Fe II (relative to those of S II and Zn II) suggest depletions slightly more severe than 
those typically found in warm, diffuse Galactic disk clouds (e.g., Savage &: Sembach 1996; Welty et 
al. 1999) and only slightly less severe than those in the strongest neutral components near v 25 
km s^^ (Table 9; component group A vs. component group DE). Analysis of the C I fine-structure 
excitation equilibrium (Jenkins Sz Shaya 1979) indicates that the thermal pressure {nuT) in that 
gas is about 3 x 10^ cm~^K — more than an order of magnitude higher than the value found for 
the strongest neutral components (Welty, Hobbs, & Morton 2002). While the increased thermal 
pressure might be suggestive of gas that has been shocked, the moderately severe depletions would 
seem to indicate that the accompanying dust grains have not been seriously disturbed. 
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5.4. Dust Survival in High- and Intermediate- Velocity Gas 

It is now well known that the depletions of various refractory elements are generally less severe 
in gas at higher LSR velocities. The initial indications of such variations in depletion came from 
observations of increased iV(Ca II)/iV(Na I) ratios at higher I^lsrI (Routly & Spitzer 1952; Sembach 
& Danks 1994). Since both Ca II and Na I are trace species in those relatively low column density 
clouds, differences in ionization could in principle produce similar effects, however. Observations 
of UV absorption lines from the dominant ions (in H I gas) of various severely depleted and mildly 
depleted elements subsequently provided strong evidence that the depletions do indeed vary with 
velocity (e.g., Shull, York, & Hobbs 1977; Jenkins et al. 1976, 1984). The relatively low spectral 
resolution and/or S/N achievable with Copernicus and lUE, and (in many cases) the incomplete 
coverage of potentially significant ions for some elements, however, generally did not allow very 
precise characterization of the depletions. 

The increase in the gas phase abundances of the refractory elements at higher IulsrI is usually 
attributed to the (partial) destruction of dust grains in interstellar shocks, primarily in the warm, 
relatively low density "intercloud" medium. A number of papers have explored the various processes 
affecting dust grains subjected to interstellar shocks — including shattering and vaporization of 
grains due to grain-grain collisions and thermal and non-thermal sputtering of grain surfaces due 
to gas-grain collisions (e.g.. Barlow k, Silk 1977; Cowie 1978; Draine k, Salpeter 1979; Scab & Shull 
1983; Jones ct al. 1994; Jones, Ticlcns, & HoUcnbach 1996). Jones et al. (1994, 1996), for example, 
find that grain destruction (i.e., the return of atoms or ions to the gas phase) peaks at a swept up 
hydrogen column density of lO^'^^-lO"'^^ cm~^, where the gas has cooled to about 10^ K (comparable 
to the values found for the HV gas toward (" Ori and 23 Ori). For a 100 km s^^ shock, the fraction 
of silicate dust destroyed (primarily by sputtering) ranges from < 10% for grains with initial radius 
< 200 A to about 30% for grains with initial radius > 1000 A. In general, the models yield time 
scales for grain destruction of order a few times 10^ years — much shorter than the time scales 
estimated for the injection of significant quantities of new grains into the ISM (e.g., in the winds 
of late- type stars). The abundances of Mg, Si, and Fe found for different interstellar environments 
(e.g., "cold, dense cloud" vs. "warm, diffuse cloud"), may suggest that the time scale for the return 
of silicon to the gas phase is significantly shorter, however (Tielens 1998). In order to maintain the 
depletions inferred for "typical" low- velocity interstellar clouds, the dust grains must therefore be 
re-assembled in the ISM — presumably via accretion of atoms and/or ions onto the grain cores that 
survived the shock(s) — during periods when the gas and dust cycle through the denser clouds. 

Determination of column densities for several significant ionization states of C, N, S, Al, Si, 
and Fe in the HV and IV components toward Q Ori affords an opportunity to gauge directly how 
much dust has survived in that ionized, perhaps shocked gas — with fewer (uncertain) corrections 
for unobserved ions required. As discussed above (§ 4.1), the HV column densities of C II + C III, 
N II -I- N III, and S II -|- S III are consistent with solar relative abundances of C, N, and S and an 
H II column density of about 7.6 x 10^^ cm~^ (though carbon may be slightly depleted). The IV 
column densities are about a factor of 3 higher (§ 4.2). The corresponding column density sums 
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for Al, Si, and Fe yield depletions of roughly —0.8 to —0.9 dex, —0.3 to —0.4 dex, and —0.95 dex, 
respectively, when compared to the inferred HV and IV iV(H II), though (as noted above) the 
depletions of aluminum and iron would be less severe (by perhaps 0.3-0.5 dex) if Al IV and Fe IV 
are significant. The residual depletions in the HV and IV gas toward ^ Ori are less severe than 
those found for warm, diffuse Galactic clouds, but may be slightly more severe than those found 
for clouds in the Galactic halo (e.g.. Savage k, Sembach 1996; Fitzpatrick 1996; Welty et al. 1999). 

If the abundances in the HV components are representative of post-shock gas, then we must 
estimate the abundances in the pre-shock gas in order to determine how much of the dust was 
destroyed in the shock. It is unlikely that the depletions in the prc-shock gas were as severe as those 
found in cold, dense clouds, given the relatively low nn and A^(H). It would not be unreasonable, 
however, to adopt the depletions characterizing warm, low-density clouds for the pre-shock gas, 
since very similar values are found for the low Af(H) components seen at low velocities toward both 
Q Ori [e.g., near km — "component 1" of Jenkins &; Peimbert (1996) or component group 
A in Table 9] and 23 Ori [the "weak, low- velocity" gas described by Welty et al. (1999)]. If those 
depletions (listed in Table 9) arc representative of the prc-shock material toward C Ori, then the 
shock appears to have returned of order 10% of the Al, Si, and Fe (and roughly 60% of the carbon) 
from the dust to the gas phase. While the ^ 10% destruction percentage for Al, Si, and Fe is 
comparable to that predicted by the theoretical models for a 100 km s~^ shock (Jones et al. 1994, 
1996), the similar values for silicon and iron stand in contrast to those derived from comparisons 
of cold cloud and warm cloud depletions that suggest that silicon is returned to the gas faster than 
iron (Tielens 1998; Jones 1999). The milder depictions of aluminum and silicon in the HV gas 
toward 23 Ori suggest that an even higher fraction of the dust was destroyed there. 

In light of the residual depletions of Al, Si, and Fe in the HV gas, the high gas-phase abundance 
of carbon (and the corresponding inferred high dust destruction percentage) is intriguing. In lower 
velocity clouds, the depletion of carbon is mild —0.4 dex, relative to solar reference abundances), 
but roughly invariant — with no discernible dependence on average sightline density < nn > or 
molecular fraction /(H2) (over the somewhat limited ranges sampled so far) (Sofia et al. 1997). 
The resulting amount of carbon in the dust is consistent with that needed to produce the various 
features in extinction and emission that have been attributed to carbonaceous material (e.g., Mathis 
1996). The depletions of most other elements (including Al, Si, and Fe) are typically more severe, 
and exhibit a clear dependence on < nn >, on /(H2), and also (evidently) on Ii'lsrI- Models for 
the destruction of dust in shocks predict that the time scale for destruction of graphite grains is 
roughly 1.5-2.0 times as long as the corresponding time scale for silicate grains (e.g., Jones et al. 
1994, 1996). We would thus expect that silicon would be returned to the gas phase more rapidly 
and more completely than carbon — but that does not appear to be the case for the HV gas toward 
C, Ori (though it may be for the HV gas toward 23 Ori). 



6. Conclusions / Summary 



We have used near-UV spectra obtained with the Hubble Space Telescope GHRS together with 
far-UV spectra from IMAPS and Copernicus to study the abundances and physical conditions in 
the ionized high- and intermediate-velocity gas seen along the line of sight to the star ( Ori A. 
The combination of high spectral rcsohition (FWHM ~ 3.3-4.5 km s^^), high S/N (90 550 for 
GHRS cchclle; 30-80 for IMAPS), and broad spectral coverage has enabled a much more detailed 
and accurate view of the HV and IV gas than was possible using spectra from Copernicus alone. 
With 2a detection limits better than 1 mA in many cases, we have obtained HV column densities 
for both singly and doubly ionized C, N, Al, Si, S, and Fe, and also for Mg II and C II*; we also 
have significantly restrictive limits for a number of neutral and triply ionized species. In the IV 
gas, wc also detect absorption from O I, C IV, and Si IV. The column density sums V(X II) + 
V(X III) for C, N, and S yield estimates for the total hydrogen column density V(H I) + iV(H II); 
by comparison, the corresponding sums for Al, Si, and Fe then provide estimates for the depletions. 
Limits (or values) for A^(0 I) and A^(N I) yield estimates for the amount of neutral gas iV(H I). Fits 
to the high-resolution line profiles suggest (at least) 5 components for the HV gas and 8 components 
for the IV gas; comparisons of line widths, for species of different atomic weight, allow estimates 
for the temperature and turbulent velocity in each component. Once T is known, analysis of the 
C II fine-structure excitation provides estimates for the local densities (ng ~ na) in several of the 
components. 

The HV gas (—105 km ^ v ^ —65 km s~^) is predominantly ionized, with 
iV(H l)/iV(H II) < 0.001 and N{X IIl)/iV(X II) ranging from 1.4 (nitrogen) to 11 (iron). The 
total hydrogen column density A^(H II) is about 7.6 x 10^^ cm^^. The depletions [—0.8 dex (Al), 
—0.4 dex (Si), and —0.95 dex (Fe), relative to C, N, and S] appear to be intermediate between 
those found in Galactic halo clouds and in warm, diffuse disk clouds — though the depletions for 
aluminum and iron would be less severe if significant amounts of Al IV and Fe IV are present. While 
the various N(X III)/V(X II) ratios would be consistent with collisional ionization equilibrium at 
temperatures of 25,000-80,000 K, the component line widths imply T ~ 9000 K — i.e., the HV gas 
is not in CIE. Characteristic densities for the HV components are rig ~ ran ~ 0.1-0.2 cm~^; typical 
thermal pressures are 2njiT ~ 2000-4000 cm' 

For the IV gas (—60 km s~^ ^ v ^ —10 km s~^), the overall column densities are all roughly 
3 times higher than for the HV gas, but the ion ratios, depletions, and inferred physical conditions 
are all similar to those found for the HV gas. 

We have considered two sets of shock models to try to understand the origins of the HV 
and IV gas toward ^ Ori: (1) for the HV gas, a shock moving toward us at a velocity of order 
-100 km s-\ and (2) for the IV shock moving away from us, produced when a flow of 

high-velocity gas encounters lower velocity material. The models arc I-D, assume steady flow, and 
include the effects of magnetic fields, stellar ionizing radiation, and radiation emitted as the gas 
behind the shock cools. The inclusion of a stellar radiation field seems necessary to reproduce the 
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observed ion ratios A'"(X III)/A'"(X II) and the weakness of N I. While the models can yield ion 
abundances and absorption-line profiles that resemble those observed, there are some significant 
differences. For example, the models have difficulty reproducing the full set of observed ion ratios; 
the broad, complex tails in the observed profiles; and the general weakness of the observed C IV 
and Si IV. Inhomogeneities in the density and/or magnetic field may produce the more complex 
structure seen in the observed profiles. More complex shock models (e.g., for slow mode MHD 
shocks and/or switch-on, switch-off shocks) may be required to more accurately reproduce the 
observed line profiles and ion ratios. 

We have also briefly examined other models for ionized gas — diffuse photoionized gas, denser 
H II regions, collisionally ionized gas in equilibrium, and isochorically cooling gas. Those diverse 
models predict very different values for the various ion ratios [e.g., A^(X III)/A^(X II)]. While 
photoionization by an 09 star radiation field might yield a reasonable match with a number of 
the HV gas ion ratios, none of the various models appears able to reproduce the full set of ratios 
observed in the HV or IV gas toward ( Ori. 

If we assume that the depletions in the pre-shock gas were similar to those found in warm, 
diffuse disk clouds — such as some of those seen at lower velocities toward ( Ori and 23 Ori — 
then the residual depletions of Al, Si, and Fe found in the HV and IV gas suggest that of order 10% 
of the dust might have been destroyed in the shocks that produced the HV and IV gas. While the 
~ 10% destruction for Al, Si, and Fe is comparable to the percentages predicted by the theoretical 
models for a 100 km shock, the similar values for silicon and iron disagree with expectations 
derived from comparisons of cold cloud and warm cloud depletions — which suggest that silicon 
is returned to the gas faster than iron. Moreover, the observed near-solar gas phase abundance of 
carbon appears inconsistent with the predicted longer time scales for the destruction of graphite 
grains. 

Many of the observations reported here were obtained with the NASA/ESA Hubble Space Tele- 
scope at the Space Telescope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS5-26555. DEW is grateful for support 
provided by NASA grant GO-05896.01-94A (administered by STScI) and by NASA LTSA grant 
NAG5-3228. Support for flying IMAPS on the ORFEUS-SPAS-I mission was provided by NASA 
grant NAG5-616 to Princeton University. The ORFEUS-SPAS project was a joint undertaking of 
the US and German space agencies NASA and DARA. 

A. High- Velocity Gas toward r CMa 

During Cycle 5, we also obtained new GHRS ECH-A spectra of the 09 lb star r CMa (V = 
4.40; E{B - F) = 0.13; d ~ 1 kpc), supplementing the ECH-B and G160M spectra described by 
Trapero et al. (1996), in order to probe in more detail the high- velocity gas along a line of sight 



-26- 



somewhat removed from the Orion-Eridanus region, at (Z, 6) = (238?2, — 5?5). Six of the ECH-A 

settings Hsted in Table 1 (centered near 1190, 1206, 1260, 1334, 1402, and 1670 A) were used, 
and tlic resulting spectra were processed similarly to those obtained for ( Ori. Since r CMa is 
significantly fainter than (" Ori, the S/Ns in the final combined, normalized spectra were lower, 
ranging from about 55 to 125. 

The normalized echelle spectra of species showing HV absorption, including the ECH-B spec- 
trum of Mg II A2796 from Trapero et al. (1996), are shown in Figure 9. In the figure (as in this 
paper) we focus on the high-velocity gas. Table 10 lists the equivalent widths or 2a upper limits for 
the high-velocity gas absorption for the transitions measured toward r CMa (including some from 
Trapero et al.). The last column in the Table gives the column density derived for each species, 
either from fits to the observed line profiles or, for upper limits, from the weak line approximation. 
There does not appear to be significant absorption at intermediate velocities toward r CMa, even 
in the strong C II A1334 and Si III A1206 lines. We note, however, that IV absorption may be 
present for C III toward r CMa, and is present for C II, C III, and Si III toward 29 CMa (24 arcmin 
away), in Copernicus spectra displayed by Cohn & York (1977). 

Table 11 lists the component parameters derived from fits to the stronger HV absorption 
features of C II, Mg II, Si II, and Si III. The three first ions are likely to be very similarly distributed, 
and so their line profiles were fitted simultaneously. The resulting velocity structure was used to 
fit the Si III line at 1206 A, though there may be slight relative shifts for several of the Si III 
components. 

While we do not have as complete a set of column densities as for ^ Ori, the significantly 
smaller A''(Si IIl)/A/'(Si II) ratio for the HV gas toward r CMa (especially in component 1) suggests 
that neither C III nor N III is likely to be significant. We may thus use the observed log[A^(C II)] 
= 13.55 cm~^ to estimate log[Ar(H II)] ~ 17.0 cm~^ — nearly an order of magnitude smaller than 
toward (, Ori and 23 Ori. The sum of the HV column densities of Si II and Si III imply that 
log[A''(H II)] > 16.96 cm~^, while the corresponding sum for Al II and Al III yields 16.64 cm~^ 
< log[A^(H II)] < 17.05 cm~^. We thus find that silicon is essentially undepleted in the HV gas 
toward r CMa (as found for the HV gas toward 23 Ori) and that aluminum is depleted by less 
than about 0.4 dex. Furthermore, A^(Mg II) (alone) suggests that log[A^(H II)] > 16.64 cm~^, so 
that magnesium is also depleted by less than about 0.4 dex and that Ar(Mg III)/Ar(Mg II) is much 
smaller in the HV gas toward r CMa than toward the two Orion stars. 

The 6-values obtained in the profile fits to the C II, Mg II, and Si II lines suggest that the 
temperatures of the individual components are similar to those in the HV gas toward ( Ori and 
23 Ori (i.e., 7000-9000 K). Since we have only an upper limit for A''(C II*), we may only derive 
corresponding limits for the local densities, pressures, and thicknesses of the HV components — 
which are generally consistent with the corresponding values found for the HV gas toward 23 Ori 
and C Ori (Table 6). 
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Fig. la. — Spectra of UV absorption lines observed toward ( Ori. The source of each profile is 
noted below the line identification. Tick marks denote the components derived from fits to the 
spectra. The vertical dotted lines separate the high- velocity (HV), intermediate- velocity (IV), and 
low-velocity (LV) component groups. Absorption due to other species is noted above some of the 
spectra; the dotted line in the N III spectrum shows the absorption due to Si II, and the dotted 
line in the S III spectrum shows the absorption due to Si II. Note the expanded vertical scale for 
some of the weaker lines; where two lines are plotted for the HV gas, the lower has been expanded 
by a factor of 5 (C II*, Al II, S II, S III, Fe II) or of 2.5 (Al III). 
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Fig. 2. — Schematic diagram of the interstellar absorption observed toward C, Ori. The asterisks 
denote the individual components found in fits to the line profiles of the various species noted 
at the right. The high- and intermediate-velocity gas (HV and IV) is predominantly ionized; the 
low- velocity gas (LV) is predominantly neutral. 
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sapadg 

Fig. 3. — Schematic diagram showing column densities of ionization states of various elements ob- 
served in the high- and intermediate-velocity gas toward ^ Ori, arranged versus ionization potential. 
Vertical dotted lines at 14.5 and 33.5 eV show the limits for species not detected in HV gas (N I 
has the highest Xion for any undetected neutral species; Si III has the lowest Xion for any species 
with undetected next higher ion); O I, C IV, and Si IV are detected in IV gas, however. 
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Fig. 4. — Ion ratios observed for HV gas toward ( Ori (filled circles; with l-cr uncertainties), 

compared with predictions from various models of ionized gas. The top panel shows the ratios 
predicted by Gruenwald & Viegas (1992) for H II regions ionized by an 04 star radiation field 
(open squares) or by a BO star field (open triangles). In each case, the upper point represents a 
sightline through the center of the region (fractional radius r = 0.0); the lower is for a sightline at 
r = 0.7; single points are for r = 0.0. The open circles show the values predicted for the "standard 
composite" model of Sembach et al. (2000) for diffuse ionized gas. The middle panel shows the 
ratios predicted by the collisional ionization equilibrium models of Sutherland & Dopita (1993), for 
T = 20,000, 40,000, and 80,000 K. The bottom panel shows the ratios predicted by two models of 
cooling (non-CIE) gas. The open squares give the ratios from Shapiro & Moore (1976); the open 
triangles give values from Schmutzler & Tscharnuter (1993) — both for T = 10^ K (the lowest T 
considered) . 
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Fig. 5. — Line profiles predicted by models of shocks moving toward us at 95-100 km s , versus 
high- velocity line profiles observed toward Ori (solid lines). The models described in Table 8 (A, 
B, C) are given by dotted, short-dashed, and long-dashed lines, respectively. Note the expanded 
vertical scale for some weak features. N I and S IV are not computed in the models, but are 
expected to be weak. The model profiles have not been convolved with the relevant instrumental 
response functions. 
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Fig. 6. — Ion ratios observed for HV gas toward Ori (filled circles; with l-a uncertainties), 
compared with predictions from several shock models (see Tables 6 and 8 and Figure 5). Model A 
is the "basic" shock model; model B adds a stronger stellar radiation field; model C has both the 
stronger radiation field and a stronger magnetic field. 
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Fig. 7. — Line profiles predicted by models of shocks moving away from us (produced by impact 
of high- velocity flow on low- velocity material), versus intermediate- velocity line profiles observed 
toward ^ Ori (solid lines). The models described in Table 8 (D, E, F) are given by dotted, short- 
dashed, and long-dashed lines, respectively. Note the expanded vertical scale for some weak features. 
While profiles for N I and S IV are not computed in the models, the total N(N I) is predicted to 
be much higher than observed for models D and E, but very low for model F; A^(S IV) is very low 
in all three models. Some of the absorption near IV S III is due to HV Si II. The model profiles 
have not been convolved with the relevant instrumental response functions. 
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Fig. 8. — Ion ratios observed for IV gas toward C Ori (filled circles; with 1-a uncertainties), 
compared with predictions from several (reverse) shock models (see Tables 6 and 8 and Figure 
7). Model D is the "basic" shock model (with no stellar radiation field); model E adds a stronger 
magnetic field; model F has both the stronger magnetic field and a strong radiation field. 
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Fig. 9. — High velocity absorption observed toward r CMa. Tick marks denote the components 
derived from fits to the spectra. Note the expanded vertical scale for Al II. 
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Table 1. C Ori — GHRS Spectra 



Order 


•^min Amax 


FWHM'i 


S/N 


Species'' 




(A) 


(km s-i) 






EchcUc Spectra 


A-47 


1188.2-1194.8 


3.56 


130-210 


Si II, Si II*, S III, C I, CI I 


A-47 


1201.2-1207.5 


3.40 


115 


Si III 


A-45 


1257.5-1264.1 


3.36 


175-265 


S II, Si II, Fc II, C I, C I*, C I** 


A-43 


1300.0-1307.2 


3.53 


235-285 


P II, O I, S I, Si II 


A-42 


1330.4-1337.8 


3.54 


220-250 


C II, C II*, CI I, P III 


A-40 


1398.3-1406.1 


3.53 


270-280 


Si IV, S I, Sn II 


A-36 


1545.3-1554.1 


3.61 


90-245 


C IV, Co II, Si I 


A-34 


1666.7-1675.5 


3.35 


105-120 


AI II, P I, Mg I 


B-31 


1801.3-1811.5 


3.56 


195-270 


Ni II, S I, Si II 


B-30 


1848.4-1859.1 


3.66 


190 


Al III, Fc I 


B-27 


2053.8-2065.7 


3.67 


315-415 


Cr II, Zn II, Co II 


B-25 


2253.3-2265.5 


3.43 


455-525 


Fe II, Al I, Ni I, Fe I 


B-24 


2334.4-2346.8 


3.50 


120 


Fe II (spybal) 


B-22 


2590.7-2604.0 


3.27 


385-555 


Mn II, Fe II 


B-20 


2790.9-2806.7 


3.55 


415-500 


Mg II, Mn I 


G160M Spectra^ 




1186.6-1224.9 


18.6 


135-225 


Si II, S III, C I, Mn II, N I, Si III, H I 




1205.7-1245.0 


18.3 


185-335 


H I, Ge II, N V, Mg II 




1220.3-1258.7 


18.0 


180-345 


Ge II, N V, Mg II, S II 




1239.5-1277.6 


17.7 


60-85 


Mg II, N V, S II, Si II, Fe II, C I, C I*, C I**, Si II* 




1315.4-1353.5 


16.5 


60-70 


C II, C II*, CI I, C I, C I*, C I**, Ni II 




1367.7-1405.7 


15.8 


35-60 


Ni II, Si IV 




1382.3-1420.1 


15.6 


35-55 


Si IV, Ga II 



^For echelle spectra, resolution = 78000 + 6.4(mA — 55000), where m is order number. For G160M 
spectra, resolution = 14400 + 16(A - 1100). 

''Possible; not all detected. 

''Prom HST archive; obtained by GO program 2584. 
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Table 2. ^ Ori — High- Velocity Gas Equivalent Widths 



Ion 


^vac 

(A) 




log(a)'' 


(mA) 




C I 


1193.0308 


E 


1.731 


<0.5 




C II 


1036.3367 


I 


2.106 


49.2±3.4 






1334.5323 


E 


2.229 


82.4±0.5 




C II* 


1335.6627 


E 


1.229 


blend with C II* 




1335.7077 


E 


2.183 


1.1±0.6 


B 


C III 


977.026 


C 


2.872 


(161.5) 




C IV 


1548.195 


E 


2.470 


<0.8 






1550.770 


E 


2.169 


<0.4 




N I 


1134.9803 


I 


1.691 


<2.0 






1199.5496 


G 


2.193 


<i.o 




N II 


1083.990 


I 


2.096 


34.6±3.3 




N III 


989.799 


C 


2.023 


(42.9) 




N V 


1238.821 


G 


2.286 


<0.5 




O I 


1302.1685 


E 


1.831 


<U.4 




Mgl 


1668.4228 


E 


0.302 


<U.O 




Mg II 


2796.352 


E 


3.234 


19.2±0.6 






2803.531 


E 


2.933 


Q O-Un T 




Al I 


2264.1637 


E 


2.305 


<0.2 




AI II 


1670.7874 


E 


3.486 


4.3±1.4 




Al III 


1854.7164 


E 


3.017 


4.3±1.0 




Si I 


1552.948 


E 


1.302 


<0.3 




Si II 


1260.4221 


E 


3.104 


24.3±0.7 






1304.3702 


E 


2.067 


2.3±0.5 




Si II* 


1264.7377 


G 


3.058 


<2.4 




Si III 


1206.500 


E 


3.304 


[95.1] 




Si IV 


1402.770 


E 


2.554 


<0.25 




S I 


1807.3113 


E 


2.239 


<0.5 




S II 


1259.519 


E 


1.311 


0.7±0.4 




S III 


1190.208 


E 


1.421 


2.1±0.7 




S IV 


1062.671 


I 


1.628 


<1.7 




Ar I 


1066.6599 


I 


1.851 


<1.4 




Mn II 


2594.499 


E 


2.847 


<0.5 




Fc 1 


2260.2102 


E 


1.170 


<0.3 




Fc II 


2600.1729 


E 


2.765 


2.9±0.4 




Fe III 


1122.526 


IC 


1.947 


3±1 





* Vacuum wavelengths and /-values from Morton (1991), 
updated by Welty et al. 1999. E = GHRS echelle; G = 
GHRS G160M; I = IMAPS; C = Copernicus. 

= major contributor to blend. Uncertainties and 
limits are 2-a. Values in parentheses are from Cowie et al. 
1979. Value for Si III in square braces is from profile fit. 
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Table 3. C Ori — HV, IV Component Structure^ 



"0 

(km 


C II 
(10^=*) 


C II* 

(10") 


(10^^) 


(10") 


(Ifll") 


('iQl") 


Si II 
(10") 


Si III 
(10^^) 


S II 
(10") 


(10^^) 


(lo") 


-97.0±0.2 
IH 


0.1B±0.04 
(3.6) 


<0.6 
(3.5) 


6.6±0.8 
(3.3) 


0.6±0.1 
(2.2) 


0.6±0.6 
(2.1) 


2.2±1.1 
(2.1) 


2.4±0.2 
(2.0) 


0.28±0.16 
(2.5) 


<3.1 
(1.9) 


<0.3 
(1.9) 


<1.9 
(1.5) 


-92.r)±().i 

2H 


3.7±0.1 
3.8±0.4 


4.4±0.5 
(3.8) 


11.4±0.6 
3.4±0.5 


2.4±0.1 
2.5±0.4 


4.2±0.4 
(2.5) 


10.1±0.9 
(2.5) 


10.8±0.2 
2.4±0.3 


3.1±0.2 
3.7±0.2 


12.3±2.8 
(2.3) 


3.4±a.3 
(2.3) 


13.2±1.9 
2.1±0.4 


-86.6±0.2 
3H 


2.0±0.1 


<0.6 


8.6±0.3 


l.ldiO.l 


3.0±0.6 


9.0±1.1 


6.4±0.2 
W-OJ 


3.0±0.1 

1^4. { ) 


7.7±3.7 

1,4. yj 


2.3±0.4 

fd Q1 
^4.y; 


4.9±2.6 

1,4. OJ 


-77.5±0.2 
4H 


1.6±0.1 
(5.5) 


0.9±0.4 
(5.5) 


8.7±0.4 
(5.5) 


0.7±0.1 
(3.9) 


1.7±0.4 
(3.7) 


3.3±0.9 
(3.7) 


2.4±0.1 
(3.4) 


2.2±0.2 
4.4±0.3 


6.3±3.0 
(3.3) 


2.2±0.3 
(3.3) 


3.7±2.0 
(2.9) 


-67.7±0.6 
5H 


0.06±0.01 
(3.5) 


<0.4 

r3.5) 


1.2±0.2 
(3.3) 


<0.06 
(2.2) 


<0.3 
(2.1) 


<0.8 
(2.1) 


<().08 
(2.0) 


0.5±0.1 
(5.0) 


<2.4 
(1.9) 


<0.3 
(1.9) 


<1.6 
(1.5) 


-68.8±0.9 
61 


0.06±0.01 

(3.5) 














l.OiO.l 

4.2±0.5 








-51.8±0.5 
71 


0.20±0.05 
2.5±0.7 














l.OiO.l 
(3.5) 








-46.9±1.0 
81 


0.36±0.04 
(3.5) 














0.8±0.2 
(3.0) 








-39.9±0.2 
91 


O.BiO.l 
3.2±0.4 














1.4±0.4 
3.1±0.7 








-32.0±0.4 
101 


1.6±0.1 

(4.2) 


<5.0 

(4.2) 


1.3±0.4 

(4.0) 


<0.04 

(2.0) 


<0.7 

(2.0) 


<1.3 

(2.0) 


1.3±1.0 

(2.0) 


2.2±0.1 

3.5±0.4 


<4.8 
(2.0) 


1.5±0.8 
(2.0) 


<2.0 

(1.9) 


-24.3±0.4 
111 


4.0±0.7 
(3.3) 


[3.9±0.3] 
(3.3) 


(5.9±0.2) 
(3.0) 


1.0±0.2 
(2.5) 


2.7±1.0 
(2.5) 


2.8±0.8 
(2.5) 


5.0±1.0 
(2.5) 


(3.1±0.2) 
(2.5) 


(2.5) 
(2.5) 


4.5±0.9 
(2.5) 


6.1±2.6 
(2.3) 


-17.8±0.4 
121 


25.5±4.4 
(3.7) 


[24.3±1.7] 
(3.7) 


(40. ±2.) 
(3.5) 


7.7±0.3 
(3.0) 


19.5±3.8 
(3.0) 


61. ±3. 
(4.0) 


32. ±2. 
(3.0) 


(18.7±1.0) 
(3.0) 


48.±7. 
(3.0) 


28. ±2. 
(3.0) 


44. ±7. 
(2.8) 


-14.1±().4 
131 


18.4±3.2 
(4.0) 


[17.5±1.2] 
(4.0) 


(30. ±1.) 
(3.8) 


4.2±0.3 
(3.5) 


13.5±4.4 
(3.5) 


16. ±3. 
(3.5) 


23. ±3. 
(3.5) 


(9.4±0.5) 
(3.5) 


15.4±7.5 
(3.5) 


6.0±2.5 
(3.5) 


27.±8. 
(3.2) 



'^Entries are N (x the factor under the ion name, in cm ^) and b (km s "'"). HV, IV values in parentheses were fixed in the fitting; values in square 
brackets were varied together. Uncertainties are l-cr; limits are 2-<t. Components (1—5) are HV; components (6—13) are IV. 
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Table 4. ^ Ori — HV, IV Column Densities and Depletions^ 



Ion 


log[N(HV)] 
(cm-2) 


log[N(IV)] 
(cm^-^) 


H I 
H II 


<14.84 


ir 01 +0.08 
io.oi_o IQ 

18 44+0.05 


C I 
C II 
C II* 
C III 

C IV 

D(C) 


<11.9 

i3.88t;j:»i 

11 70+0. 07 

14 1+°-^ 
<11.47 
-n lo+'^-i^ 

'^■^^-0.13 


<11.9 

14.711;;:°^ 

12 66+0-02 

14 7+0.2 
' -0.2 

+0.0210-0 


N I 
N II 

N III 
N V 
D(N) 


<12.0 
10 cr+O.Ol 

10.00_Q QJ 

1 ■? 7+0.1 
' -0.1 

<11.4 

+0.081°:°^ 


<12.0 
13.89troi 

<11.4 
+0.03l°:fo 


O I 


<11.71 


i2.i8i°;fo 


Mgl 
Mg II 


<13.3 

- - ——1-0 09 

ll.D( _g (J2 


<13.3 
19 11 ^^"^.^-^ 


Al I 
Al II 
Al III 
D(A1) 


<10.7 
10.98int 

-0.82l°;°i 


<10.7 
11.55i«:°? 


Si I 
Si II 
Si II* 
Si III 
Si IV 
D(Si) 


<12.0 

ir, qn + O.Ol 

<11.3 

12.951^:; 
<10.93 


<12.0 
12 7Q+0 02 

<11.3 

,0 c:7+0.08 
1,5.0 (_(, „() 

-0.34l°:fo 


S I 
S II 
S III 

S IV 
D(S) 


<11.3 

12 47+0 09 
^2 Q1+0.03 

<12.65 
-0 11+004 


<11.3 

19 S9 + 0.06 
13 60+0 03 

<12.65 
-0 04+0-05 


Fe I 
Fe II 

Fc III 
D(Fe) 


<12.0 

12.38+0 17 

U.y _o.l6 


<12.0 

11 8S+0.05 



^Unccrtainitcs arc l-cr; limits are 2- 
cr. H I estimated from O I; H II es- 
timated from C II + C III, N II + 
N III, S II + S III — using reference 
abundances from Anders & Grevesse 
1989 and Grevesse & Noels 1993. De- 
pletions are relative to those reference 
abundances. 
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Table 5. ^ Ori — HV, IV Component Properties 



Comp* N{C n*)/N{C 11)^ Ue = n-a" log[Af(H 11)]^ r log{2nHT) 

(cm~3) (cm~2) (pc) (cm~3K) 



IH 


<0.04 


<0.7 


16.5 


>0.015 


<4.3 


2H 


0.012±0.001 


0.21±0.03 


17.48 


0.47 


3.6 


3H 


<0.003 


<0.05 


17.36 


>1.5 


<3.0 


4H 


0.006±0.003 


0.10±0.05 


17.25 


0.58 


3.3 


5H 


<0.07 


<1.2 


16.3 


>0.005 


<4.3 


61 


<0.08 


<1.5 


16.7 


>0.01 


<4.4 


71 


<0.025 


<0.44 


16.75 


>0.04 


<3.9 


81 


<0.014 


<0.25 


16.65 


>0.06 


<3.7 


91 


<0.006 


<0.11 


16.9 


>0.2 


<3.3 


101 


<0.03 


<0.55 


17.1 


>0.07 


<4.0 


111 


[0.010±0.002] 


[0.17±0.04] 


[16.95] 


[0.17] 


[3.5] 


121 


[0.010±0.002] 


[0.17±0.04] 


[18.15] 


[2.7] 


[3.5] 


131 


[0.010±0.002] 


[0.17±0.04] 


[17.75] 


[1.1] 


[3.5] 



Components (1-5) are HV; components (6-13) are IV. 

''Ratio assumed equal for components 111—131 in fits to C II* profile. 

"= Assuming T = 9000±2000 K, from 6-values for C II, N II, Mg II, Si II, and Fe II 
in component 2H. Limits for components 61—91 assume iV(C II*) < 5 X lO'^" cm~^. 

dProm iV(C II) + A'(C III), Ar(N II) + Ar(N III), Ar(S II) + Ar(S III), and/or /V(Si 
II) + N{S\ III) assuming solar reference abundances, adjusted for adopted total HV 
or IV log[iV(H II)] (see Table 4). 
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Table 6. HV and IV Gas toward C Ori, 23 Ori, r CMa^ 









HV Models 






IV Models 


Property 


C Ori HV 


23 Ori HV'^ r CMa HV 




C Ori IV 


23 Ori IV'' 


D/E/pd 



log[Ar(H II)] (cm-2) 




17.8±0.1 


IT.OitO.l 


N(n i)/N(n II) 


<0.001 


<0.01 


<0.05 


N{C III)/W{C II) 


1.7±0.8 


<0.4 




n\c IV)/N{C 11) 


<0.004 


<0.004 


<0.019 


iV(N HI)/Af(N II) 


1.4±0.3 


<0.1 




iV(Al III)/Af(Al II) 


2.6±0.3 


O.SitO.l 


<1.5 


Ar(Si iii)/Ar(Si n) 


4.4±1.0 


3.0±0.5 


0.7±0.1 


AT (Si IV)/iV(Si II) 


<0.04 


<0.06 


<0.16 


iV(S III)/iV(S II) 


2.7±0.6 


<0.6 




iV(Fe III)/iV(Fe II) 


11. ±4. 






D(C) 


-O.litO.l 


>-0.2 


(0.0) 


D(A1) 


-0.8±0.05 


-0.5±0.1 


>-0.4 


D(Si) 


-0.4±0.1 


-O.litO.l 


0.0±0.05 


D(Fe) 


-0.95±0.15 






log[Ar(C ll*)/N{C II)] 




1 co+0.05 
^•°'^-0.04 


<-1.95 


< ne > (cm'^)"" 


n 14+0.05 

'-'■.'-^-0.04 




<0.18 


log[2< riH >T] (cm-^K) 


3.4±0.2 


3.9±0.1 


<3.5 


Ar(H II)/< nn > (pc) 


1.8 


0.4 


>0.2 





18.44±0.05 


17.85±0.15 






0.001 


<0.02 




0.7/6.6/3.6 


1.0±0.5 




0.12/0.23/4.9 


0.02/0.11/0.0005 


0.011±0.002 




0.011/0.013/0.17 


1.0/2.0/0.5 


2.6±0.7 




0.16/0.23/1.8 




2.2±0.4 


1.0±0.7 




17/240/141 


6.2±1.3 


4.4±2.1 


0.13/0.34/537 


2.2/34/10 


0.13±0.01 




0.031/0.054/49 


1.6/42/20 


6.2±1.1 




0.44/1.0/91 


52/117/71 










O.OitO.l 








-0.85±0.05 


-0.35±0.15 






-0.35±0.1 


-0.15±0.15 





16+°°* 1 l^+0-44 

3.45±0.2 4.25±0.2 
5.1 0.2 



^Total or mean values over all high- or intermediate-velocity components. Uncertainties are l-u; limits are 2-<t. Depletions are relative to 
reference abundances of Anders & Grevesse 1989 and Grevesse & Noels 1993. 
''Welty et al. 1999. 

'^Shock travelling toward us at 95-100 km s^^ (§5.3.1; Tabic 7). 

''Shock travelling away from us, produced by impact of 95-100 km s^-' flow upon low-velocity material (§5.3.2; Table 7). 
<= Assumes T = 9000±2000 K for C Ori, T = 8000±1000 K for 23 Ori and r CMa. 
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Table 7. C Ori — HV Ion Ratios^ 



Element (X/S)© [X II/S II] [X III/S III] [X/S] 



c 


1.28 


+0.131°;?? 


-0.09l°;i^ 


-o.02l°;l^ 


N 


0.70 


+0.381°;?? 


+0.09l°;?g 


+0.i9l°;°f 


Al 


-0.79 


-0.70l°;?? 


^•'^-0.05 


r. 79+0.06 
'J- '^-0.05 


Si 


0.28 


_0 44+0-09 


r, 90+0.09 

'^■^'^-0.12 


-0.28l°:?« 


Fc 


n.24 


-1 07+0.11 

-0,14 


"■ ' ' -0.18 


-0.87li-« 



'^Uncertainties are l-ir; columns 3-5 are relative to solar ratios 
given in column 2 (Anders & Grevesse 1989; Grevesse & Noels 
1993). 
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Table 8. Shock Model Parameters 



Model 




Vq no -Bo 






(km s-i) 


(km s-i) (cm-3) (^G) 


(erg cm~^ A~^) 


lligli-V(>lo(:ity Gas 


A 


95 


-12 0.0032 0.5 


1 X lO"'^ 


B 


95 


-14 0.0032 0.5 


1 X 10-5 


C 


100 


-38 0.01 2.5 


1 X 10-5 


Intermediate- Velocity Gas 


D 


95 


-14 0.010 0.2 





E 


100 


-7 0.020 1.1 





F 


100 


-7 0.020 1.1 


5.1 X 10^5 



"^For diluted Kurucz (1979) model spectrum, near 1000 A. 
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Table 9. Depletions and Dust Destruction 



Element 




Galactic'^ 
Warm/Cold 


23 Ori'' 
HV/WLV/SLV 


C Ori= 
HV/A/DE 


% Destroyed'^ 


C 


8.55 


-0.4/-0.4 


> -0.2/-0.31/-0.31 


-0.12/ ■■■ / ••• 


60-67 


N 


7.97 


-0.1/-0.1 


••• /-0.15/-0.15 


+0.08/ ••■ / ••• 




Al 


6.48 


-1.1/-2.4 


-0.5/-1.36/-1.71 


-0.82/-1.43/-1.84 


9-13 


Si 


7.55 


-0.4/-1.3 


-0.1/-0.52/-0.89 


-0.39/-0.59/-0.72 


0-20 


Fe 


7.51 


-1.4/-2.2 


••• /-1.49/-1.92 


-0.97/-1.46/-1.85 


6-8 



^Logarithmic solcir abundances from Anders & Grevesse 1989 and Grevesse & Noels 1993. H = 
12.0. 

''Welty et al. 1999; WLV = "weak, low-velocity" components with log[Af(H)] = 19.6 cm-^; SLV 
= "strong, low-velocity" components with log[Ar(H)] = 20.7 cm~^. 

"^Welty et al., in prep.; A = components near km s~^; DE = strong components near +25 km 



Derived from f Ori H V (post-shock) versus C Ori A or Galactic warm cloud (assumed pre-shock) 
values. 
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Table 10. r CMa — HV Gas Equivalent Widths and Column Densities 



Ion 


\ 3- 

^vac 

(A) 




log(fA)^ 


(mA) 


log{N) 
(cm- 2) 


C I 


1560.3092 


G 


2.050 


<1.9° 


<12.09 


C II 
C II* 


1334.5323 

1335.6627 


E 

E 


2.229 

1.229 


53.6±1.2 

blend with C 11* 


io cc+O.Ol 
lO.O0_Q 02 

<11.60 




iooo. / U / / 


T? 

tj 


z. loo 


<U.o rS 




C IV 


1548.195 


G 


2.470 


<2.8= 


<11.84 




1550.770 


G 


2.169 


<2.8'= 




N V 


1238.821 


G 


2.286 


<3.7 


<12.23 




1242.804 


G 


1.986 


<3.5 




O I 


1302.1685 


G 


1.831 


<3.2= 


<12.61 


A/To- T 


ZUZU.^ / OO 


E 


z.ooy 


^9 7^ 


<^ i-i-.Oi- 


Mg II 


2796.352 


E 


3.234 


58.6±4.6'= 


19 99+0.04 




2803.531 


E 


2.933 


25.1±4.4<= 




A 1 T 


2601 . ( lOZ 


Hi 


z.oyo 


<l.o 


<ii.04 


Al II 


1670.7874 


E 


3.486 


5.9±2.4 


11.11_0.22 


Al III 


1854.7164 


E 


3.017 


<3.4'= 


<11.30 




1862.7895 


E 


2.716 


<3.4>= 




Si I 


2515.0725 


E 


2.723 


<2.6'^ 


<11.34 


Si II 


1190.4158 
1193.2897 


E 
E 


2.474 
2.775 


5.5±1.7 
10.9±1.7 


12 so+o ol 




1260.4221 


E 


3.104 


23.3±1.5 




Si II* 


1264.7377 


G 


3.058 


<2.4 


<11.28 


Si III 
Si IV 


1206.500 
1393.755 


E 
G 


3.304 
2.855 


24.9±2.3 
<2.7= 


19 19+0.03 
<11.49 




1402.770 


E 


2.554 


<2.7'= 




S I 


1807.3113 


E 


2.239 


<1.7'= 


<11.79 


S II 


1259.519 


E 


1.311 


<0.9 


<12.60 


S III 


1190.208 


E 


1.421 


<1.4 


<12.70 


Fe I 


2484.0211 


E 


3.141 


<2.0= 


<10.82 


Fe II 


2374.4612 


E 


1.871 


<2.7'= 


<12.24 



^Vacuum wavelengths and /-values from Morton (1991), updated by 
Welty ct al. 1999. E = GHRS echelle; G = G160M. 

'^B = major contributor to blend. Uncertainties and limits are 2-a. 

■^Values from Trapero at al. 1996. 
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Tablc 11. T CMa HV Component Stnicture'" 



Vq 


C II 


Mg II 


Si II 


Si III 


(km s-l) 


(1012) 


(lOii) 


(10") 


(lOii) 


-96.4±0.1 


12.6±1.1 


7.0±0.3 


9.5±0.4 


1.3±0.3 


IH 


3.3±0.2 


2.3±0.2 


2.4±0.2 


1.8±0.9 


-87.8±0.1 


17.3±1.6 


7.3±0.4 


8.5±0.6 


8.4±0.6 


2H 


5.2±0.6 


3.8±0.3 


4.1±0.4 


5.0±0.6 


-80.9±0.4 


1.9±0.7 


0.9±0.2 


(0.2) 


1.2±0.4 


3H 


(3.1) 


(2.1) 


(2.0) 


(2.5) 


-69.7±0.2 


4.2±0.2 


1.5±0.2 


1.8±0.2 


2.4±0.3 


4H 


4.9±0.4 


2.9±0.5 


2.4±0.5 


4.3±0.7 



''Entries are N (cm ^) and b (km s i). Values in paren- 
theses were fixed in the fitting. Uncertainties are 1-a. 



